      module mo_lu_factor
      private
      public :: lu_fac
      contains
      subroutine lu_fac01( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(1) = 1._r8 / lu(1)
         lu(2) = 1._r8 / lu(2)
         lu(3) = 1._r8 / lu(3)
         lu(4) = 1._r8 / lu(4)
         lu(5) = 1._r8 / lu(5)
         lu(6) = lu(6) * lu(5)
         lu(7) = 1._r8 / lu(7)
         lu(8) = 1._r8 / lu(8)
         lu(9) = lu(9) * lu(8)
         lu(10) = 1._r8 / lu(10)
         lu(11) = 1._r8 / lu(11)
         lu(12) = 1._r8 / lu(12)
         lu(13) = 1._r8 / lu(13)
         lu(14) = 1._r8 / lu(14)
         lu(15) = 1._r8 / lu(15)
         lu(16) = 1._r8 / lu(16)
         lu(17) = 1._r8 / lu(17)
         lu(18) = 1._r8 / lu(18)
         lu(19) = 1._r8 / lu(19)
         lu(20) = lu(20) * lu(19)
         lu(634) = lu(634) - lu(20) * lu(567)
         lu(21) = 1._r8 / lu(21)
         lu(22) = lu(22) * lu(21)
         lu(23) = lu(23) * lu(21)
         lu(616) = lu(616) - lu(22) * lu(568)
         lu(634) = lu(634) - lu(23) * lu(568)
         lu(24) = 1._r8 / lu(24)
         lu(25) = lu(25) * lu(24)
         lu(26) = lu(26) * lu(24)
         lu(617) = lu(617) - lu(25) * lu(569)
         lu(634) = lu(634) - lu(26) * lu(569)
         lu(27) = 1._r8 / lu(27)
         lu(28) = lu(28) * lu(27)
         lu(29) = lu(29) * lu(27)
         lu(582) = lu(582) - lu(28) * lu(570)
         lu(634) = lu(634) - lu(29) * lu(570)
         lu(30) = 1._r8 / lu(30)
         lu(31) = lu(31) * lu(30)
         lu(32) = lu(32) * lu(30)
         lu(615) = lu(615) - lu(31) * lu(571)
         lu(634) = lu(634) - lu(32) * lu(571)
         lu(34) = 1._r8 / lu(34)
         lu(35) = lu(35) * lu(34)
         lu(77) = lu(77) - lu(35) * lu(74)
         lu(634) = lu(634) - lu(35) * lu(572)
         lu(796) = lu(796) - lu(35) * lu(770)
         lu(36) = 1._r8 / lu(36)
         lu(37) = lu(37) * lu(36)
         lu(38) = lu(38) * lu(36)
         lu(39) = lu(39) * lu(36)
         lu(40) = lu(40) * lu(36)
         lu(574) = lu(574) - lu(37) * lu(573)
         lu(601) = lu(601) - lu(38) * lu(573)
         lu(634) = lu(634) - lu(39) * lu(573)
         lu(635) = lu(635) - lu(40) * lu(573)
         lu(41) = 1._r8 / lu(41)
         lu(42) = lu(42) * lu(41)
         lu(43) = lu(43) * lu(41)
         lu(578) = lu(578) - lu(42) * lu(574)
         lu(634) = lu(634) - lu(43) * lu(574)
         lu(44) = 1._r8 / lu(44)
         lu(45) = lu(45) * lu(44)
         lu(46) = lu(46) * lu(44)
         lu(47) = lu(47) * lu(44)
         lu(48) = lu(48) * lu(44)
         lu(586) = lu(586) - lu(45) * lu(575)
         lu(600) = lu(600) - lu(46) * lu(575)
         lu(634) = lu(634) - lu(47) * lu(575)
         lu(635) = lu(635) - lu(48) * lu(575)
         lu(49) = 1._r8 / lu(49)
         lu(50) = lu(50) * lu(49)
         lu(51) = lu(51) * lu(49)
         lu(634) = lu(634) - lu(50) * lu(576)
         lu(635) = lu(635) - lu(51) * lu(576)
         lu(684) = lu(684) - lu(50) * lu(641)
         lu(685) = lu(685) - lu(51) * lu(641)
         lu(52) = 1._r8 / lu(52)
         lu(53) = lu(53) * lu(52)
         lu(54) = lu(54) * lu(52)
         lu(55) = lu(55) * lu(52)
         lu(141) = - lu(53) * lu(139)
         lu(142) = - lu(54) * lu(139)
         lu(144) = - lu(55) * lu(139)
         lu(537) = lu(537) - lu(53) * lu(526)
         lu(552) = lu(552) - lu(54) * lu(526)
         lu(560) = lu(560) - lu(55) * lu(526)
      end subroutine lu_fac01
      subroutine lu_fac02( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(56) = 1._r8 / lu(56)
         lu(57) = lu(57) * lu(56)
         lu(58) = lu(58) * lu(56)
         lu(451) = lu(451) - lu(57) * lu(444)
         lu(453) = - lu(58) * lu(444)
         lu(628) = lu(628) - lu(57) * lu(577)
         lu(634) = lu(634) - lu(58) * lu(577)
         lu(678) = lu(678) - lu(57) * lu(642)
         lu(684) = lu(684) - lu(58) * lu(642)
         lu(59) = 1._r8 / lu(59)
         lu(60) = lu(60) * lu(59)
         lu(61) = lu(61) * lu(59)
         lu(62) = lu(62) * lu(59)
         lu(602) = - lu(60) * lu(578)
         lu(635) = lu(635) - lu(61) * lu(578)
         lu(638) = lu(638) - lu(62) * lu(578)
         lu(753) = lu(753) - lu(60) * lu(750)
         lu(764) = lu(764) - lu(61) * lu(750)
         lu(767) = lu(767) - lu(62) * lu(750)
         lu(63) = 1._r8 / lu(63)
         lu(64) = lu(64) * lu(63)
         lu(65) = lu(65) * lu(63)
         lu(66) = lu(66) * lu(63)
         lu(67) = lu(67) * lu(63)
         lu(68) = lu(68) * lu(63)
         lu(69) = lu(69) * lu(63)
         lu(696) = lu(696) - lu(64) * lu(691)
         lu(707) = lu(707) - lu(65) * lu(691)
         lu(713) = lu(713) - lu(66) * lu(691)
         lu(714) = lu(714) - lu(67) * lu(691)
         lu(715) = lu(715) - lu(68) * lu(691)
         lu(717) = lu(717) - lu(69) * lu(691)
         lu(70) = 1._r8 / lu(70)
         lu(71) = lu(71) * lu(70)
         lu(72) = lu(72) * lu(70)
         lu(73) = lu(73) * lu(70)
         lu(755) = lu(755) - lu(71) * lu(751)
         lu(767) = lu(767) - lu(72) * lu(751)
         lu(768) = lu(768) - lu(73) * lu(751)
         lu(776) = lu(776) - lu(71) * lu(771)
         lu(800) = lu(800) - lu(72) * lu(771)
         lu(801) = lu(801) - lu(73) * lu(771)
         lu(75) = 1._r8 / lu(75)
         lu(76) = lu(76) * lu(75)
         lu(77) = lu(77) * lu(75)
         lu(78) = lu(78) * lu(75)
         lu(79) = lu(79) * lu(75)
         lu(604) = lu(604) - lu(76) * lu(579)
         lu(634) = lu(634) - lu(77) * lu(579)
         lu(635) = lu(635) - lu(78) * lu(579)
         lu(639) = lu(639) - lu(79) * lu(579)
         lu(776) = lu(776) - lu(76) * lu(772)
         lu(796) = lu(796) - lu(77) * lu(772)
         lu(797) = lu(797) - lu(78) * lu(772)
         lu(801) = lu(801) - lu(79) * lu(772)
         lu(80) = 1._r8 / lu(80)
         lu(81) = lu(81) * lu(80)
         lu(82) = lu(82) * lu(80)
         lu(83) = lu(83) * lu(80)
         lu(322) = lu(322) - lu(81) * lu(318)
         lu(325) = - lu(82) * lu(318)
         lu(326) = lu(326) - lu(83) * lu(318)
         lu(621) = lu(621) - lu(81) * lu(580)
         lu(634) = lu(634) - lu(82) * lu(580)
         lu(635) = lu(635) - lu(83) * lu(580)
         lu(731) = lu(731) - lu(81) * lu(721)
         lu(743) = - lu(82) * lu(721)
         lu(744) = lu(744) - lu(83) * lu(721)
         lu(84) = 1._r8 / lu(84)
         lu(85) = lu(85) * lu(84)
         lu(86) = lu(86) * lu(84)
         lu(87) = lu(87) * lu(84)
         lu(88) = lu(88) * lu(84)
         lu(404) = lu(404) - lu(85) * lu(402)
         lu(407) = lu(407) - lu(86) * lu(402)
         lu(409) = - lu(87) * lu(402)
         lu(413) = lu(413) - lu(88) * lu(402)
         lu(623) = lu(623) - lu(85) * lu(581)
         lu(627) = lu(627) - lu(86) * lu(581)
         lu(634) = lu(634) - lu(87) * lu(581)
         lu(640) = lu(640) - lu(88) * lu(581)
         lu(673) = lu(673) - lu(85) * lu(643)
         lu(677) = lu(677) - lu(86) * lu(643)
         lu(684) = lu(684) - lu(87) * lu(643)
         lu(690) = lu(690) - lu(88) * lu(643)
         lu(89) = 1._r8 / lu(89)
         lu(90) = lu(90) * lu(89)
         lu(91) = lu(91) * lu(89)
         lu(92) = lu(92) * lu(89)
         lu(93) = lu(93) * lu(89)
         lu(94) = lu(94) * lu(89)
         lu(95) = lu(95) * lu(89)
         lu(540) = lu(540) - lu(90) * lu(527)
         lu(546) = lu(546) - lu(91) * lu(527)
         lu(552) = lu(552) - lu(92) * lu(527)
         lu(558) = lu(558) - lu(93) * lu(527)
         lu(560) = lu(560) - lu(94) * lu(527)
         lu(563) = lu(563) - lu(95) * lu(527)
         lu(614) = lu(614) - lu(90) * lu(582)
         lu(621) = lu(621) - lu(91) * lu(582)
         lu(627) = lu(627) - lu(92) * lu(582)
         lu(633) = - lu(93) * lu(582)
         lu(635) = lu(635) - lu(94) * lu(582)
         lu(638) = lu(638) - lu(95) * lu(582)
         lu(96) = 1._r8 / lu(96)
         lu(97) = lu(97) * lu(96)
         lu(98) = lu(98) * lu(96)
         lu(99) = lu(99) * lu(96)
         lu(100) = lu(100) * lu(96)
         lu(271) = lu(271) - lu(97) * lu(270)
         lu(272) = lu(272) - lu(98) * lu(270)
         lu(274) = - lu(99) * lu(270)
         lu(277) = lu(277) - lu(100) * lu(270)
         lu(612) = lu(612) - lu(97) * lu(583)
         lu(621) = lu(621) - lu(98) * lu(583)
         lu(634) = lu(634) - lu(99) * lu(583)
         lu(640) = lu(640) - lu(100) * lu(583)
         lu(663) = lu(663) - lu(97) * lu(644)
         lu(671) = - lu(98) * lu(644)
         lu(684) = lu(684) - lu(99) * lu(644)
         lu(690) = lu(690) - lu(100) * lu(644)
      end subroutine lu_fac02
      subroutine lu_fac03( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(101) = 1._r8 / lu(101)
         lu(102) = lu(102) * lu(101)
         lu(103) = lu(103) * lu(101)
         lu(104) = lu(104) * lu(101)
         lu(105) = lu(105) * lu(101)
         lu(482) = lu(482) - lu(102) * lu(474)
         lu(483) = - lu(103) * lu(474)
         lu(485) = - lu(104) * lu(474)
         lu(486) = lu(486) - lu(105) * lu(474)
         lu(630) = lu(630) - lu(102) * lu(584)
         lu(632) = lu(632) - lu(103) * lu(584)
         lu(634) = lu(634) - lu(104) * lu(584)
         lu(635) = lu(635) - lu(105) * lu(584)
         lu(680) = lu(680) - lu(102) * lu(645)
         lu(682) = lu(682) - lu(103) * lu(645)
         lu(684) = lu(684) - lu(104) * lu(645)
         lu(685) = lu(685) - lu(105) * lu(645)
         lu(106) = 1._r8 / lu(106)
         lu(107) = lu(107) * lu(106)
         lu(108) = lu(108) * lu(106)
         lu(109) = lu(109) * lu(106)
         lu(110) = lu(110) * lu(106)
         lu(111) = lu(111) * lu(106)
         lu(112) = lu(112) * lu(106)
         lu(586) = lu(586) - lu(107) * lu(585)
         lu(592) = lu(592) - lu(108) * lu(585)
         lu(627) = lu(627) - lu(109) * lu(585)
         lu(634) = lu(634) - lu(110) * lu(585)
         lu(635) = lu(635) - lu(111) * lu(585)
         lu(636) = lu(636) - lu(112) * lu(585)
         lu(693) = lu(693) - lu(107) * lu(692)
         lu(694) = - lu(108) * lu(692)
         lu(707) = lu(707) - lu(109) * lu(692)
         lu(714) = lu(714) - lu(110) * lu(692)
         lu(715) = lu(715) - lu(111) * lu(692)
         lu(716) = lu(716) - lu(112) * lu(692)
         lu(113) = 1._r8 / lu(113)
         lu(114) = lu(114) * lu(113)
         lu(115) = lu(115) * lu(113)
         lu(160) = - lu(114) * lu(156)
         lu(161) = lu(161) - lu(115) * lu(156)
         lu(559) = lu(559) - lu(114) * lu(528)
         lu(560) = lu(560) - lu(115) * lu(528)
         lu(634) = lu(634) - lu(114) * lu(586)
         lu(635) = lu(635) - lu(115) * lu(586)
         lu(684) = lu(684) - lu(114) * lu(646)
         lu(685) = lu(685) - lu(115) * lu(646)
         lu(714) = lu(714) - lu(114) * lu(693)
         lu(715) = lu(715) - lu(115) * lu(693)
         lu(116) = 1._r8 / lu(116)
         lu(117) = lu(117) * lu(116)
         lu(118) = lu(118) * lu(116)
         lu(119) = lu(119) * lu(116)
         lu(120) = lu(120) * lu(116)
         lu(330) = lu(330) - lu(117) * lu(329)
         lu(331) = lu(331) - lu(118) * lu(329)
         lu(336) = - lu(119) * lu(329)
         lu(337) = lu(337) - lu(120) * lu(329)
         lu(614) = lu(614) - lu(117) * lu(587)
         lu(617) = lu(617) - lu(118) * lu(587)
         lu(634) = lu(634) - lu(119) * lu(587)
         lu(635) = lu(635) - lu(120) * lu(587)
         lu(665) = - lu(117) * lu(647)
         lu(669) = lu(669) - lu(118) * lu(647)
         lu(684) = lu(684) - lu(119) * lu(647)
         lu(685) = lu(685) - lu(120) * lu(647)
         lu(121) = 1._r8 / lu(121)
         lu(122) = lu(122) * lu(121)
         lu(123) = lu(123) * lu(121)
         lu(124) = lu(124) * lu(121)
         lu(125) = lu(125) * lu(121)
         lu(634) = lu(634) - lu(122) * lu(588)
         lu(635) = lu(635) - lu(123) * lu(588)
         lu(638) = lu(638) - lu(124) * lu(588)
         lu(639) = lu(639) - lu(125) * lu(588)
         lu(684) = lu(684) - lu(122) * lu(648)
         lu(685) = lu(685) - lu(123) * lu(648)
         lu(688) = lu(688) - lu(124) * lu(648)
         lu(689) = lu(689) - lu(125) * lu(648)
         lu(763) = lu(763) - lu(122) * lu(752)
         lu(764) = lu(764) - lu(123) * lu(752)
         lu(767) = lu(767) - lu(124) * lu(752)
         lu(768) = lu(768) - lu(125) * lu(752)
         lu(126) = 1._r8 / lu(126)
         lu(127) = lu(127) * lu(126)
         lu(128) = lu(128) * lu(126)
         lu(129) = lu(129) * lu(126)
         lu(130) = lu(130) * lu(126)
         lu(321) = lu(321) - lu(127) * lu(319)
         lu(322) = lu(322) - lu(128) * lu(319)
         lu(325) = lu(325) - lu(129) * lu(319)
         lu(326) = lu(326) - lu(130) * lu(319)
         lu(616) = lu(616) - lu(127) * lu(589)
         lu(621) = lu(621) - lu(128) * lu(589)
         lu(634) = lu(634) - lu(129) * lu(589)
         lu(635) = lu(635) - lu(130) * lu(589)
         lu(668) = lu(668) - lu(127) * lu(649)
         lu(671) = lu(671) - lu(128) * lu(649)
         lu(684) = lu(684) - lu(129) * lu(649)
         lu(685) = lu(685) - lu(130) * lu(649)
         lu(131) = 1._r8 / lu(131)
         lu(132) = lu(132) * lu(131)
         lu(133) = lu(133) * lu(131)
         lu(134) = lu(134) * lu(131)
         lu(135) = lu(135) * lu(131)
         lu(627) = lu(627) - lu(132) * lu(590)
         lu(634) = lu(634) - lu(133) * lu(590)
         lu(635) = lu(635) - lu(134) * lu(590)
         lu(637) = lu(637) - lu(135) * lu(590)
         lu(677) = lu(677) - lu(132) * lu(650)
         lu(684) = lu(684) - lu(133) * lu(650)
         lu(685) = lu(685) - lu(134) * lu(650)
         lu(687) = lu(687) - lu(135) * lu(650)
         lu(736) = lu(736) - lu(132) * lu(722)
         lu(743) = lu(743) - lu(133) * lu(722)
         lu(744) = lu(744) - lu(134) * lu(722)
         lu(746) = lu(746) - lu(135) * lu(722)
         lu(136) = 1._r8 / lu(136)
         lu(137) = lu(137) * lu(136)
         lu(138) = lu(138) * lu(136)
         lu(344) = - lu(137) * lu(341)
         lu(345) = lu(345) - lu(138) * lu(341)
         lu(498) = - lu(137) * lu(491)
         lu(504) = - lu(138) * lu(491)
         lu(553) = lu(553) - lu(137) * lu(529)
         lu(559) = lu(559) - lu(138) * lu(529)
         lu(628) = lu(628) - lu(137) * lu(591)
         lu(634) = lu(634) - lu(138) * lu(591)
         lu(737) = lu(737) - lu(137) * lu(723)
         lu(743) = lu(743) - lu(138) * lu(723)
         lu(790) = lu(790) - lu(137) * lu(773)
         lu(796) = lu(796) - lu(138) * lu(773)
         lu(812) = lu(812) - lu(137) * lu(803)
         lu(818) = - lu(138) * lu(803)
         lu(140) = 1._r8 / lu(140)
         lu(141) = lu(141) * lu(140)
         lu(142) = lu(142) * lu(140)
         lu(143) = lu(143) * lu(140)
         lu(144) = lu(144) * lu(140)
         lu(145) = lu(145) * lu(140)
         lu(537) = lu(537) - lu(141) * lu(530)
         lu(552) = lu(552) - lu(142) * lu(530)
         lu(558) = lu(558) - lu(143) * lu(530)
         lu(560) = lu(560) - lu(144) * lu(530)
         lu(563) = lu(563) - lu(145) * lu(530)
         lu(609) = lu(609) - lu(141) * lu(592)
         lu(627) = lu(627) - lu(142) * lu(592)
         lu(633) = lu(633) - lu(143) * lu(592)
         lu(635) = lu(635) - lu(144) * lu(592)
         lu(638) = lu(638) - lu(145) * lu(592)
         lu(697) = - lu(141) * lu(694)
         lu(707) = lu(707) - lu(142) * lu(694)
         lu(713) = lu(713) - lu(143) * lu(694)
         lu(715) = lu(715) - lu(144) * lu(694)
         lu(718) = lu(718) - lu(145) * lu(694)
      end subroutine lu_fac03
      subroutine lu_fac04( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(146) = 1._r8 / lu(146)
         lu(147) = lu(147) * lu(146)
         lu(148) = lu(148) * lu(146)
         lu(149) = lu(149) * lu(146)
         lu(150) = lu(150) * lu(146)
         lu(151) = lu(151) * lu(146)
         lu(194) = lu(194) - lu(147) * lu(193)
         lu(195) = lu(195) - lu(148) * lu(193)
         lu(196) = lu(196) - lu(149) * lu(193)
         lu(197) = lu(197) - lu(150) * lu(193)
         lu(199) = - lu(151) * lu(193)
         lu(600) = lu(600) - lu(147) * lu(593)
         lu(601) = lu(601) - lu(148) * lu(593)
         lu(602) = lu(602) - lu(149) * lu(593)
         lu(626) = lu(626) - lu(150) * lu(593)
         lu(634) = lu(634) - lu(151) * lu(593)
         lu(657) = - lu(147) * lu(651)
         lu(658) = lu(658) - lu(148) * lu(651)
         lu(659) = - lu(149) * lu(651)
         lu(676) = - lu(150) * lu(651)
         lu(684) = lu(684) - lu(151) * lu(651)
         lu(152) = 1._r8 / lu(152)
         lu(153) = lu(153) * lu(152)
         lu(154) = lu(154) * lu(152)
         lu(155) = lu(155) * lu(152)
         lu(310) = - lu(153) * lu(301)
         lu(313) = - lu(154) * lu(301)
         lu(316) = lu(316) - lu(155) * lu(301)
         lu(373) = - lu(153) * lu(367)
         lu(376) = lu(376) - lu(154) * lu(367)
         lu(380) = - lu(155) * lu(367)
         lu(551) = lu(551) - lu(153) * lu(531)
         lu(559) = lu(559) - lu(154) * lu(531)
         lu(563) = lu(563) - lu(155) * lu(531)
         lu(626) = lu(626) - lu(153) * lu(594)
         lu(634) = lu(634) - lu(154) * lu(594)
         lu(638) = lu(638) - lu(155) * lu(594)
         lu(788) = lu(788) - lu(153) * lu(774)
         lu(796) = lu(796) - lu(154) * lu(774)
         lu(800) = lu(800) - lu(155) * lu(774)
         lu(157) = 1._r8 / lu(157)
         lu(158) = lu(158) * lu(157)
         lu(159) = lu(159) * lu(157)
         lu(160) = lu(160) * lu(157)
         lu(161) = lu(161) * lu(157)
         lu(162) = lu(162) * lu(157)
         lu(438) = lu(438) - lu(158) * lu(436)
         lu(439) = - lu(159) * lu(436)
         lu(440) = lu(440) - lu(160) * lu(436)
         lu(441) = lu(441) - lu(161) * lu(436)
         lu(442) = - lu(162) * lu(436)
         lu(552) = lu(552) - lu(158) * lu(532)
         lu(558) = lu(558) - lu(159) * lu(532)
         lu(559) = lu(559) - lu(160) * lu(532)
         lu(560) = lu(560) - lu(161) * lu(532)
         lu(563) = lu(563) - lu(162) * lu(532)
         lu(677) = lu(677) - lu(158) * lu(652)
         lu(683) = lu(683) - lu(159) * lu(652)
         lu(684) = lu(684) - lu(160) * lu(652)
         lu(685) = lu(685) - lu(161) * lu(652)
         lu(688) = lu(688) - lu(162) * lu(652)
         lu(163) = 1._r8 / lu(163)
         lu(164) = lu(164) * lu(163)
         lu(165) = lu(165) * lu(163)
         lu(166) = lu(166) * lu(163)
         lu(167) = lu(167) * lu(163)
         lu(516) = lu(516) - lu(164) * lu(511)
         lu(519) = - lu(165) * lu(511)
         lu(522) = lu(522) - lu(166) * lu(511)
         lu(525) = lu(525) - lu(167) * lu(511)
         lu(627) = lu(627) - lu(164) * lu(595)
         lu(634) = lu(634) - lu(165) * lu(595)
         lu(637) = lu(637) - lu(166) * lu(595)
         lu(640) = lu(640) - lu(167) * lu(595)
         lu(677) = lu(677) - lu(164) * lu(653)
         lu(684) = lu(684) - lu(165) * lu(653)
         lu(687) = lu(687) - lu(166) * lu(653)
         lu(690) = lu(690) - lu(167) * lu(653)
         lu(811) = lu(811) - lu(164) * lu(804)
         lu(818) = lu(818) - lu(165) * lu(804)
         lu(821) = lu(821) - lu(166) * lu(804)
         lu(824) = lu(824) - lu(167) * lu(804)
         lu(168) = 1._r8 / lu(168)
         lu(169) = lu(169) * lu(168)
         lu(170) = lu(170) * lu(168)
         lu(361) = lu(361) - lu(169) * lu(349)
         lu(364) = - lu(170) * lu(349)
         lu(376) = lu(376) - lu(169) * lu(368)
         lu(379) = lu(379) - lu(170) * lu(368)
         lu(519) = lu(519) - lu(169) * lu(512)
         lu(522) = lu(522) - lu(170) * lu(512)
         lu(634) = lu(634) - lu(169) * lu(596)
         lu(637) = lu(637) - lu(170) * lu(596)
         lu(684) = lu(684) - lu(169) * lu(654)
         lu(687) = lu(687) - lu(170) * lu(654)
         lu(714) = lu(714) - lu(169) * lu(695)
         lu(717) = lu(717) - lu(170) * lu(695)
         lu(743) = lu(743) - lu(169) * lu(724)
         lu(746) = lu(746) - lu(170) * lu(724)
         lu(818) = lu(818) - lu(169) * lu(805)
         lu(821) = lu(821) - lu(170) * lu(805)
         lu(171) = 1._r8 / lu(171)
         lu(172) = lu(172) * lu(171)
         lu(173) = lu(173) * lu(171)
         lu(174) = lu(174) * lu(171)
         lu(175) = lu(175) * lu(171)
         lu(261) = - lu(172) * lu(259)
         lu(264) = - lu(173) * lu(259)
         lu(267) = lu(267) - lu(174) * lu(259)
         lu(269) = - lu(175) * lu(259)
         lu(304) = - lu(172) * lu(302)
         lu(307) = - lu(173) * lu(302)
         lu(313) = lu(313) - lu(174) * lu(302)
         lu(317) = - lu(175) * lu(302)
         lu(538) = lu(538) - lu(172) * lu(533)
         lu(543) = lu(543) - lu(173) * lu(533)
         lu(559) = lu(559) - lu(174) * lu(533)
         lu(565) = lu(565) - lu(175) * lu(533)
         lu(612) = lu(612) - lu(172) * lu(597)
         lu(616) = lu(616) - lu(173) * lu(597)
         lu(634) = lu(634) - lu(174) * lu(597)
         lu(640) = lu(640) - lu(175) * lu(597)
         lu(176) = 1._r8 / lu(176)
         lu(177) = lu(177) * lu(176)
         lu(178) = lu(178) * lu(176)
         lu(179) = lu(179) * lu(176)
         lu(180) = lu(180) * lu(176)
         lu(181) = lu(181) * lu(176)
         lu(182) = lu(182) * lu(176)
         lu(279) = lu(279) - lu(177) * lu(278)
         lu(280) = lu(280) - lu(178) * lu(278)
         lu(281) = - lu(179) * lu(278)
         lu(282) = lu(282) - lu(180) * lu(278)
         lu(284) = - lu(181) * lu(278)
         lu(285) = lu(285) - lu(182) * lu(278)
         lu(613) = lu(613) - lu(177) * lu(598)
         lu(621) = lu(621) - lu(178) * lu(598)
         lu(624) = lu(624) - lu(179) * lu(598)
         lu(627) = lu(627) - lu(180) * lu(598)
         lu(634) = lu(634) - lu(181) * lu(598)
         lu(635) = lu(635) - lu(182) * lu(598)
         lu(664) = lu(664) - lu(177) * lu(655)
         lu(671) = lu(671) - lu(178) * lu(655)
         lu(674) = - lu(179) * lu(655)
         lu(677) = lu(677) - lu(180) * lu(655)
         lu(684) = lu(684) - lu(181) * lu(655)
         lu(685) = lu(685) - lu(182) * lu(655)
      end subroutine lu_fac04
      subroutine lu_fac05( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(183) = 1._r8 / lu(183)
         lu(184) = lu(184) * lu(183)
         lu(185) = lu(185) * lu(183)
         lu(186) = lu(186) * lu(183)
         lu(187) = lu(187) * lu(183)
         lu(188) = lu(188) * lu(183)
         lu(189) = lu(189) * lu(183)
         lu(391) = lu(391) - lu(184) * lu(390)
         lu(392) = lu(392) - lu(185) * lu(390)
         lu(394) = lu(394) - lu(186) * lu(390)
         lu(395) = lu(395) - lu(187) * lu(390)
         lu(397) = - lu(188) * lu(390)
         lu(398) = lu(398) - lu(189) * lu(390)
         lu(614) = lu(614) - lu(184) * lu(599)
         lu(622) = lu(622) - lu(185) * lu(599)
         lu(625) = lu(625) - lu(186) * lu(599)
         lu(629) = lu(629) - lu(187) * lu(599)
         lu(634) = lu(634) - lu(188) * lu(599)
         lu(635) = lu(635) - lu(189) * lu(599)
         lu(665) = lu(665) - lu(184) * lu(656)
         lu(672) = lu(672) - lu(185) * lu(656)
         lu(675) = lu(675) - lu(186) * lu(656)
         lu(679) = lu(679) - lu(187) * lu(656)
         lu(684) = lu(684) - lu(188) * lu(656)
         lu(685) = lu(685) - lu(189) * lu(656)
         lu(190) = 1._r8 / lu(190)
         lu(191) = lu(191) * lu(190)
         lu(192) = lu(192) * lu(190)
         lu(199) = lu(199) - lu(191) * lu(194)
         lu(200) = lu(200) - lu(192) * lu(194)
         lu(205) = - lu(191) * lu(202)
         lu(206) = lu(206) - lu(192) * lu(202)
         lu(248) = lu(248) - lu(191) * lu(245)
         lu(249) = lu(249) - lu(192) * lu(245)
         lu(453) = lu(453) - lu(191) * lu(445)
         lu(454) = lu(454) - lu(192) * lu(445)
         lu(559) = lu(559) - lu(191) * lu(534)
         lu(560) = lu(560) - lu(192) * lu(534)
         lu(634) = lu(634) - lu(191) * lu(600)
         lu(635) = lu(635) - lu(192) * lu(600)
         lu(684) = lu(684) - lu(191) * lu(657)
         lu(685) = lu(685) - lu(192) * lu(657)
         lu(743) = lu(743) - lu(191) * lu(725)
         lu(744) = lu(744) - lu(192) * lu(725)
         lu(796) = lu(796) - lu(191) * lu(775)
         lu(797) = lu(797) - lu(192) * lu(775)
         lu(195) = 1._r8 / lu(195)
         lu(196) = lu(196) * lu(195)
         lu(197) = lu(197) * lu(195)
         lu(198) = lu(198) * lu(195)
         lu(199) = lu(199) * lu(195)
         lu(200) = lu(200) * lu(195)
         lu(201) = lu(201) * lu(195)
         lu(536) = lu(536) - lu(196) * lu(535)
         lu(551) = lu(551) - lu(197) * lu(535)
         lu(558) = lu(558) - lu(198) * lu(535)
         lu(559) = lu(559) - lu(199) * lu(535)
         lu(560) = lu(560) - lu(200) * lu(535)
         lu(563) = lu(563) - lu(201) * lu(535)
         lu(602) = lu(602) - lu(196) * lu(601)
         lu(626) = lu(626) - lu(197) * lu(601)
         lu(633) = lu(633) - lu(198) * lu(601)
         lu(634) = lu(634) - lu(199) * lu(601)
         lu(635) = lu(635) - lu(200) * lu(601)
         lu(638) = lu(638) - lu(201) * lu(601)
         lu(659) = lu(659) - lu(196) * lu(658)
         lu(676) = lu(676) - lu(197) * lu(658)
         lu(683) = lu(683) - lu(198) * lu(658)
         lu(684) = lu(684) - lu(199) * lu(658)
         lu(685) = lu(685) - lu(200) * lu(658)
         lu(688) = lu(688) - lu(201) * lu(658)
         lu(203) = 1._r8 / lu(203)
         lu(204) = lu(204) * lu(203)
         lu(205) = lu(205) * lu(203)
         lu(206) = lu(206) * lu(203)
         lu(207) = lu(207) * lu(203)
         lu(551) = lu(551) - lu(204) * lu(536)
         lu(559) = lu(559) - lu(205) * lu(536)
         lu(560) = lu(560) - lu(206) * lu(536)
         lu(565) = lu(565) - lu(207) * lu(536)
         lu(626) = lu(626) - lu(204) * lu(602)
         lu(634) = lu(634) - lu(205) * lu(602)
         lu(635) = lu(635) - lu(206) * lu(602)
         lu(640) = lu(640) - lu(207) * lu(602)
         lu(676) = lu(676) - lu(204) * lu(659)
         lu(684) = lu(684) - lu(205) * lu(659)
         lu(685) = lu(685) - lu(206) * lu(659)
         lu(690) = lu(690) - lu(207) * lu(659)
         lu(759) = - lu(204) * lu(753)
         lu(763) = lu(763) - lu(205) * lu(753)
         lu(764) = lu(764) - lu(206) * lu(753)
         lu(769) = lu(769) - lu(207) * lu(753)
         lu(208) = 1._r8 / lu(208)
         lu(209) = lu(209) * lu(208)
         lu(210) = lu(210) * lu(208)
         lu(211) = lu(211) * lu(208)
         lu(212) = lu(212) * lu(208)
         lu(213) = lu(213) * lu(208)
         lu(214) = lu(214) * lu(208)
         lu(627) = lu(627) - lu(209) * lu(603)
         lu(634) = lu(634) - lu(210) * lu(603)
         lu(637) = lu(637) - lu(211) * lu(603)
         lu(638) = lu(638) - lu(212) * lu(603)
         lu(639) = lu(639) - lu(213) * lu(603)
         lu(640) = lu(640) - lu(214) * lu(603)
         lu(760) = - lu(209) * lu(754)
         lu(763) = lu(763) - lu(210) * lu(754)
         lu(766) = - lu(211) * lu(754)
         lu(767) = lu(767) - lu(212) * lu(754)
         lu(768) = lu(768) - lu(213) * lu(754)
         lu(769) = lu(769) - lu(214) * lu(754)
         lu(811) = lu(811) - lu(209) * lu(806)
         lu(818) = lu(818) - lu(210) * lu(806)
         lu(821) = lu(821) - lu(211) * lu(806)
         lu(822) = lu(822) - lu(212) * lu(806)
         lu(823) = - lu(213) * lu(806)
         lu(824) = lu(824) - lu(214) * lu(806)
         lu(215) = 1._r8 / lu(215)
         lu(216) = lu(216) * lu(215)
         lu(217) = lu(217) * lu(215)
         lu(218) = lu(218) * lu(215)
         lu(384) = lu(384) - lu(216) * lu(382)
         lu(387) = - lu(217) * lu(382)
         lu(388) = lu(388) - lu(218) * lu(382)
         lu(431) = lu(431) - lu(216) * lu(429)
         lu(433) = - lu(217) * lu(429)
         lu(434) = lu(434) - lu(218) * lu(429)
         lu(440) = lu(440) - lu(216) * lu(437)
         lu(442) = lu(442) - lu(217) * lu(437)
         lu(443) = lu(443) - lu(218) * lu(437)
         lu(634) = lu(634) - lu(216) * lu(604)
         lu(638) = lu(638) - lu(217) * lu(604)
         lu(639) = lu(639) - lu(218) * lu(604)
         lu(763) = lu(763) - lu(216) * lu(755)
         lu(767) = lu(767) - lu(217) * lu(755)
         lu(768) = lu(768) - lu(218) * lu(755)
         lu(796) = lu(796) - lu(216) * lu(776)
         lu(800) = lu(800) - lu(217) * lu(776)
         lu(801) = lu(801) - lu(218) * lu(776)
         lu(219) = 1._r8 / lu(219)
         lu(220) = lu(220) * lu(219)
         lu(221) = lu(221) * lu(219)
         lu(222) = lu(222) * lu(219)
         lu(223) = lu(223) * lu(219)
         lu(224) = lu(224) * lu(219)
         lu(633) = lu(633) - lu(220) * lu(605)
         lu(634) = lu(634) - lu(221) * lu(605)
         lu(635) = lu(635) - lu(222) * lu(605)
         lu(636) = lu(636) - lu(223) * lu(605)
         lu(638) = lu(638) - lu(224) * lu(605)
         lu(683) = lu(683) - lu(220) * lu(660)
         lu(684) = lu(684) - lu(221) * lu(660)
         lu(685) = lu(685) - lu(222) * lu(660)
         lu(686) = lu(686) - lu(223) * lu(660)
         lu(688) = lu(688) - lu(224) * lu(660)
         lu(713) = lu(713) - lu(220) * lu(696)
         lu(714) = lu(714) - lu(221) * lu(696)
         lu(715) = lu(715) - lu(222) * lu(696)
         lu(716) = lu(716) - lu(223) * lu(696)
         lu(718) = lu(718) - lu(224) * lu(696)
         lu(762) = lu(762) - lu(220) * lu(756)
         lu(763) = lu(763) - lu(221) * lu(756)
         lu(764) = lu(764) - lu(222) * lu(756)
         lu(765) = lu(765) - lu(223) * lu(756)
         lu(767) = lu(767) - lu(224) * lu(756)
         lu(225) = 1._r8 / lu(225)
         lu(226) = lu(226) * lu(225)
         lu(227) = lu(227) * lu(225)
         lu(228) = lu(228) * lu(225)
         lu(229) = lu(229) * lu(225)
         lu(230) = lu(230) * lu(225)
         lu(231) = lu(231) * lu(225)
         lu(232) = lu(232) * lu(225)
         lu(514) = - lu(226) * lu(513)
         lu(516) = lu(516) - lu(227) * lu(513)
         lu(517) = lu(517) - lu(228) * lu(513)
         lu(519) = lu(519) - lu(229) * lu(513)
         lu(520) = lu(520) - lu(230) * lu(513)
         lu(523) = lu(523) - lu(231) * lu(513)
         lu(524) = lu(524) - lu(232) * lu(513)
         lu(624) = lu(624) - lu(226) * lu(606)
         lu(627) = lu(627) - lu(227) * lu(606)
         lu(632) = lu(632) - lu(228) * lu(606)
         lu(634) = lu(634) - lu(229) * lu(606)
         lu(635) = lu(635) - lu(230) * lu(606)
         lu(638) = lu(638) - lu(231) * lu(606)
         lu(639) = lu(639) - lu(232) * lu(606)
         lu(758) = - lu(226) * lu(757)
         lu(760) = lu(760) - lu(227) * lu(757)
         lu(761) = lu(761) - lu(228) * lu(757)
         lu(763) = lu(763) - lu(229) * lu(757)
         lu(764) = lu(764) - lu(230) * lu(757)
         lu(767) = lu(767) - lu(231) * lu(757)
         lu(768) = lu(768) - lu(232) * lu(757)
      end subroutine lu_fac05
      subroutine lu_fac06( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(233) = 1._r8 / lu(233)
         lu(234) = lu(234) * lu(233)
         lu(235) = lu(235) * lu(233)
         lu(236) = lu(236) * lu(233)
         lu(237) = lu(237) * lu(233)
         lu(238) = lu(238) * lu(233)
         lu(239) = lu(239) * lu(233)
         lu(240) = lu(240) * lu(233)
         lu(495) = lu(495) - lu(234) * lu(492)
         lu(497) = lu(497) - lu(235) * lu(492)
         lu(498) = lu(498) - lu(236) * lu(492)
         lu(499) = lu(499) - lu(237) * lu(492)
         lu(501) = lu(501) - lu(238) * lu(492)
         lu(504) = lu(504) - lu(239) * lu(492)
         lu(505) = lu(505) - lu(240) * lu(492)
         lu(625) = lu(625) - lu(234) * lu(607)
         lu(627) = lu(627) - lu(235) * lu(607)
         lu(628) = lu(628) - lu(236) * lu(607)
         lu(629) = lu(629) - lu(237) * lu(607)
         lu(631) = lu(631) - lu(238) * lu(607)
         lu(634) = lu(634) - lu(239) * lu(607)
         lu(635) = lu(635) - lu(240) * lu(607)
         lu(675) = lu(675) - lu(234) * lu(661)
         lu(677) = lu(677) - lu(235) * lu(661)
         lu(678) = lu(678) - lu(236) * lu(661)
         lu(679) = lu(679) - lu(237) * lu(661)
         lu(681) = lu(681) - lu(238) * lu(661)
         lu(684) = lu(684) - lu(239) * lu(661)
         lu(685) = lu(685) - lu(240) * lu(661)
         lu(241) = 1._r8 / lu(241)
         lu(242) = lu(242) * lu(241)
         lu(243) = lu(243) * lu(241)
         lu(244) = lu(244) * lu(241)
         lu(323) = lu(323) - lu(242) * lu(320)
         lu(325) = lu(325) - lu(243) * lu(320)
         lu(326) = lu(326) - lu(244) * lu(320)
         lu(407) = lu(407) - lu(242) * lu(403)
         lu(409) = lu(409) - lu(243) * lu(403)
         lu(410) = lu(410) - lu(244) * lu(403)
         lu(450) = lu(450) - lu(242) * lu(446)
         lu(453) = lu(453) - lu(243) * lu(446)
         lu(454) = lu(454) - lu(244) * lu(446)
         lu(480) = lu(480) - lu(242) * lu(475)
         lu(485) = lu(485) - lu(243) * lu(475)
         lu(486) = lu(486) - lu(244) * lu(475)
         lu(497) = lu(497) - lu(242) * lu(493)
         lu(504) = lu(504) - lu(243) * lu(493)
         lu(505) = lu(505) - lu(244) * lu(493)
         lu(627) = lu(627) - lu(242) * lu(608)
         lu(634) = lu(634) - lu(243) * lu(608)
         lu(635) = lu(635) - lu(244) * lu(608)
         lu(736) = lu(736) - lu(242) * lu(726)
         lu(743) = lu(743) - lu(243) * lu(726)
         lu(744) = lu(744) - lu(244) * lu(726)
         lu(246) = 1._r8 / lu(246)
         lu(247) = lu(247) * lu(246)
         lu(248) = lu(248) * lu(246)
         lu(249) = lu(249) * lu(246)
         lu(450) = lu(450) - lu(247) * lu(447)
         lu(453) = lu(453) - lu(248) * lu(447)
         lu(454) = lu(454) - lu(249) * lu(447)
         lu(480) = lu(480) - lu(247) * lu(476)
         lu(485) = lu(485) - lu(248) * lu(476)
         lu(486) = lu(486) - lu(249) * lu(476)
         lu(552) = lu(552) - lu(247) * lu(537)
         lu(559) = lu(559) - lu(248) * lu(537)
         lu(560) = lu(560) - lu(249) * lu(537)
         lu(627) = lu(627) - lu(247) * lu(609)
         lu(634) = lu(634) - lu(248) * lu(609)
         lu(635) = lu(635) - lu(249) * lu(609)
         lu(707) = lu(707) - lu(247) * lu(697)
         lu(714) = lu(714) - lu(248) * lu(697)
         lu(715) = lu(715) - lu(249) * lu(697)
         lu(736) = lu(736) - lu(247) * lu(727)
         lu(743) = lu(743) - lu(248) * lu(727)
         lu(744) = lu(744) - lu(249) * lu(727)
         lu(789) = lu(789) - lu(247) * lu(777)
         lu(796) = lu(796) - lu(248) * lu(777)
         lu(797) = lu(797) - lu(249) * lu(777)
         lu(811) = lu(811) - lu(247) * lu(807)
         lu(818) = lu(818) - lu(248) * lu(807)
         lu(819) = lu(819) - lu(249) * lu(807)
         lu(250) = 1._r8 / lu(250)
         lu(251) = lu(251) * lu(250)
         lu(252) = lu(252) * lu(250)
         lu(253) = lu(253) * lu(250)
         lu(254) = lu(254) * lu(250)
         lu(255) = lu(255) * lu(250)
         lu(256) = lu(256) * lu(250)
         lu(257) = lu(257) * lu(250)
         lu(258) = lu(258) * lu(250)
         lu(622) = lu(622) - lu(251) * lu(610)
         lu(625) = lu(625) - lu(252) * lu(610)
         lu(629) = lu(629) - lu(253) * lu(610)
         lu(634) = lu(634) - lu(254) * lu(610)
         lu(635) = lu(635) - lu(255) * lu(610)
         lu(636) = lu(636) - lu(256) * lu(610)
         lu(638) = lu(638) - lu(257) * lu(610)
         lu(639) = lu(639) - lu(258) * lu(610)
         lu(702) = - lu(251) * lu(698)
         lu(705) = lu(705) - lu(252) * lu(698)
         lu(709) = lu(709) - lu(253) * lu(698)
         lu(714) = lu(714) - lu(254) * lu(698)
         lu(715) = lu(715) - lu(255) * lu(698)
         lu(716) = lu(716) - lu(256) * lu(698)
         lu(718) = lu(718) - lu(257) * lu(698)
         lu(719) = lu(719) - lu(258) * lu(698)
         lu(784) = lu(784) - lu(251) * lu(778)
         lu(787) = lu(787) - lu(252) * lu(778)
         lu(791) = lu(791) - lu(253) * lu(778)
         lu(796) = lu(796) - lu(254) * lu(778)
         lu(797) = lu(797) - lu(255) * lu(778)
         lu(798) = lu(798) - lu(256) * lu(778)
         lu(800) = lu(800) - lu(257) * lu(778)
         lu(801) = lu(801) - lu(258) * lu(778)
         lu(260) = 1._r8 / lu(260)
         lu(261) = lu(261) * lu(260)
         lu(262) = lu(262) * lu(260)
         lu(263) = lu(263) * lu(260)
         lu(264) = lu(264) * lu(260)
         lu(265) = lu(265) * lu(260)
         lu(266) = lu(266) * lu(260)
         lu(267) = lu(267) * lu(260)
         lu(268) = lu(268) * lu(260)
         lu(269) = lu(269) * lu(260)
         lu(304) = lu(304) - lu(261) * lu(303)
         lu(305) = lu(305) - lu(262) * lu(303)
         lu(306) = lu(306) - lu(263) * lu(303)
         lu(307) = lu(307) - lu(264) * lu(303)
         lu(308) = lu(308) - lu(265) * lu(303)
         lu(311) = lu(311) - lu(266) * lu(303)
         lu(313) = lu(313) - lu(267) * lu(303)
         lu(314) = lu(314) - lu(268) * lu(303)
         lu(317) = lu(317) - lu(269) * lu(303)
         lu(612) = lu(612) - lu(261) * lu(611)
         lu(614) = lu(614) - lu(262) * lu(611)
         lu(615) = lu(615) - lu(263) * lu(611)
         lu(616) = lu(616) - lu(264) * lu(611)
         lu(621) = lu(621) - lu(265) * lu(611)
         lu(627) = lu(627) - lu(266) * lu(611)
         lu(634) = lu(634) - lu(267) * lu(611)
         lu(635) = lu(635) - lu(268) * lu(611)
         lu(640) = lu(640) - lu(269) * lu(611)
         lu(663) = lu(663) - lu(261) * lu(662)
         lu(665) = lu(665) - lu(262) * lu(662)
         lu(667) = lu(667) - lu(263) * lu(662)
         lu(668) = lu(668) - lu(264) * lu(662)
         lu(671) = lu(671) - lu(265) * lu(662)
         lu(677) = lu(677) - lu(266) * lu(662)
         lu(684) = lu(684) - lu(267) * lu(662)
         lu(685) = lu(685) - lu(268) * lu(662)
         lu(690) = lu(690) - lu(269) * lu(662)
         lu(271) = 1._r8 / lu(271)
         lu(272) = lu(272) * lu(271)
         lu(273) = lu(273) * lu(271)
         lu(274) = lu(274) * lu(271)
         lu(275) = lu(275) * lu(271)
         lu(276) = lu(276) * lu(271)
         lu(277) = lu(277) * lu(271)
         lu(308) = lu(308) - lu(272) * lu(304)
         lu(312) = lu(312) - lu(273) * lu(304)
         lu(313) = lu(313) - lu(274) * lu(304)
         lu(314) = lu(314) - lu(275) * lu(304)
         lu(316) = lu(316) - lu(276) * lu(304)
         lu(317) = lu(317) - lu(277) * lu(304)
         lu(546) = lu(546) - lu(272) * lu(538)
         lu(558) = lu(558) - lu(273) * lu(538)
         lu(559) = lu(559) - lu(274) * lu(538)
         lu(560) = lu(560) - lu(275) * lu(538)
         lu(563) = lu(563) - lu(276) * lu(538)
         lu(565) = lu(565) - lu(277) * lu(538)
         lu(621) = lu(621) - lu(272) * lu(612)
         lu(633) = lu(633) - lu(273) * lu(612)
         lu(634) = lu(634) - lu(274) * lu(612)
         lu(635) = lu(635) - lu(275) * lu(612)
         lu(638) = lu(638) - lu(276) * lu(612)
         lu(640) = lu(640) - lu(277) * lu(612)
         lu(671) = lu(671) - lu(272) * lu(663)
         lu(683) = lu(683) - lu(273) * lu(663)
         lu(684) = lu(684) - lu(274) * lu(663)
         lu(685) = lu(685) - lu(275) * lu(663)
         lu(688) = lu(688) - lu(276) * lu(663)
         lu(690) = lu(690) - lu(277) * lu(663)
         lu(279) = 1._r8 / lu(279)
         lu(280) = lu(280) * lu(279)
         lu(281) = lu(281) * lu(279)
         lu(282) = lu(282) * lu(279)
         lu(283) = lu(283) * lu(279)
         lu(284) = lu(284) * lu(279)
         lu(285) = lu(285) * lu(279)
         lu(286) = lu(286) * lu(279)
         lu(371) = lu(371) - lu(280) * lu(369)
         lu(372) = - lu(281) * lu(369)
         lu(374) = lu(374) - lu(282) * lu(369)
         lu(375) = - lu(283) * lu(369)
         lu(376) = lu(376) - lu(284) * lu(369)
         lu(377) = lu(377) - lu(285) * lu(369)
         lu(380) = lu(380) - lu(286) * lu(369)
         lu(546) = lu(546) - lu(280) * lu(539)
         lu(549) = lu(549) - lu(281) * lu(539)
         lu(552) = lu(552) - lu(282) * lu(539)
         lu(558) = lu(558) - lu(283) * lu(539)
         lu(559) = lu(559) - lu(284) * lu(539)
         lu(560) = lu(560) - lu(285) * lu(539)
         lu(563) = lu(563) - lu(286) * lu(539)
         lu(621) = lu(621) - lu(280) * lu(613)
         lu(624) = lu(624) - lu(281) * lu(613)
         lu(627) = lu(627) - lu(282) * lu(613)
         lu(633) = lu(633) - lu(283) * lu(613)
         lu(634) = lu(634) - lu(284) * lu(613)
         lu(635) = lu(635) - lu(285) * lu(613)
         lu(638) = lu(638) - lu(286) * lu(613)
         lu(671) = lu(671) - lu(280) * lu(664)
         lu(674) = lu(674) - lu(281) * lu(664)
         lu(677) = lu(677) - lu(282) * lu(664)
         lu(683) = lu(683) - lu(283) * lu(664)
         lu(684) = lu(684) - lu(284) * lu(664)
         lu(685) = lu(685) - lu(285) * lu(664)
         lu(688) = lu(688) - lu(286) * lu(664)
      end subroutine lu_fac06
      subroutine lu_fac07( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(287) = 1._r8 / lu(287)
         lu(288) = lu(288) * lu(287)
         lu(289) = lu(289) * lu(287)
         lu(290) = lu(290) * lu(287)
         lu(291) = lu(291) * lu(287)
         lu(309) = - lu(288) * lu(305)
         lu(313) = lu(313) - lu(289) * lu(305)
         lu(315) = - lu(290) * lu(305)
         lu(317) = lu(317) - lu(291) * lu(305)
         lu(333) = - lu(288) * lu(330)
         lu(336) = lu(336) - lu(289) * lu(330)
         lu(338) = lu(338) - lu(290) * lu(330)
         lu(340) = - lu(291) * lu(330)
         lu(393) = - lu(288) * lu(391)
         lu(397) = lu(397) - lu(289) * lu(391)
         lu(399) = - lu(290) * lu(391)
         lu(401) = - lu(291) * lu(391)
         lu(548) = lu(548) - lu(288) * lu(540)
         lu(559) = lu(559) - lu(289) * lu(540)
         lu(562) = lu(562) - lu(290) * lu(540)
         lu(565) = lu(565) - lu(291) * lu(540)
         lu(623) = lu(623) - lu(288) * lu(614)
         lu(634) = lu(634) - lu(289) * lu(614)
         lu(637) = lu(637) - lu(290) * lu(614)
         lu(640) = lu(640) - lu(291) * lu(614)
         lu(673) = lu(673) - lu(288) * lu(665)
         lu(684) = lu(684) - lu(289) * lu(665)
         lu(687) = lu(687) - lu(290) * lu(665)
         lu(690) = lu(690) - lu(291) * lu(665)
         lu(732) = lu(732) - lu(288) * lu(728)
         lu(743) = lu(743) - lu(289) * lu(728)
         lu(746) = lu(746) - lu(290) * lu(728)
         lu(749) = lu(749) - lu(291) * lu(728)
         lu(292) = 1._r8 / lu(292)
         lu(293) = lu(293) * lu(292)
         lu(294) = lu(294) * lu(292)
         lu(295) = lu(295) * lu(292)
         lu(296) = lu(296) * lu(292)
         lu(297) = lu(297) * lu(292)
         lu(298) = lu(298) * lu(292)
         lu(299) = lu(299) * lu(292)
         lu(300) = lu(300) * lu(292)
         lu(351) = - lu(293) * lu(350)
         lu(354) = lu(354) - lu(294) * lu(350)
         lu(355) = lu(355) - lu(295) * lu(350)
         lu(357) = lu(357) - lu(296) * lu(350)
         lu(360) = - lu(297) * lu(350)
         lu(362) = lu(362) - lu(298) * lu(350)
         lu(365) = - lu(299) * lu(350)
         lu(366) = lu(366) - lu(300) * lu(350)
         lu(545) = lu(545) - lu(293) * lu(541)
         lu(550) = lu(550) - lu(294) * lu(541)
         lu(552) = lu(552) - lu(295) * lu(541)
         lu(554) = lu(554) - lu(296) * lu(541)
         lu(558) = lu(558) - lu(297) * lu(541)
         lu(560) = lu(560) - lu(298) * lu(541)
         lu(563) = lu(563) - lu(299) * lu(541)
         lu(564) = lu(564) - lu(300) * lu(541)
         lu(670) = lu(670) - lu(293) * lu(666)
         lu(675) = lu(675) - lu(294) * lu(666)
         lu(677) = lu(677) - lu(295) * lu(666)
         lu(679) = lu(679) - lu(296) * lu(666)
         lu(683) = lu(683) - lu(297) * lu(666)
         lu(685) = lu(685) - lu(298) * lu(666)
         lu(688) = lu(688) - lu(299) * lu(666)
         lu(689) = lu(689) - lu(300) * lu(666)
         lu(780) = lu(780) - lu(293) * lu(779)
         lu(787) = lu(787) - lu(294) * lu(779)
         lu(789) = lu(789) - lu(295) * lu(779)
         lu(791) = lu(791) - lu(296) * lu(779)
         lu(795) = lu(795) - lu(297) * lu(779)
         lu(797) = lu(797) - lu(298) * lu(779)
         lu(800) = lu(800) - lu(299) * lu(779)
         lu(801) = lu(801) - lu(300) * lu(779)
         lu(306) = 1._r8 / lu(306)
         lu(307) = lu(307) * lu(306)
         lu(308) = lu(308) * lu(306)
         lu(309) = lu(309) * lu(306)
         lu(310) = lu(310) * lu(306)
         lu(311) = lu(311) * lu(306)
         lu(312) = lu(312) * lu(306)
         lu(313) = lu(313) * lu(306)
         lu(314) = lu(314) * lu(306)
         lu(315) = lu(315) * lu(306)
         lu(316) = lu(316) * lu(306)
         lu(317) = lu(317) * lu(306)
         lu(543) = lu(543) - lu(307) * lu(542)
         lu(546) = lu(546) - lu(308) * lu(542)
         lu(548) = lu(548) - lu(309) * lu(542)
         lu(551) = lu(551) - lu(310) * lu(542)
         lu(552) = lu(552) - lu(311) * lu(542)
         lu(558) = lu(558) - lu(312) * lu(542)
         lu(559) = lu(559) - lu(313) * lu(542)
         lu(560) = lu(560) - lu(314) * lu(542)
         lu(562) = lu(562) - lu(315) * lu(542)
         lu(563) = lu(563) - lu(316) * lu(542)
         lu(565) = lu(565) - lu(317) * lu(542)
         lu(616) = lu(616) - lu(307) * lu(615)
         lu(621) = lu(621) - lu(308) * lu(615)
         lu(623) = lu(623) - lu(309) * lu(615)
         lu(626) = lu(626) - lu(310) * lu(615)
         lu(627) = lu(627) - lu(311) * lu(615)
         lu(633) = lu(633) - lu(312) * lu(615)
         lu(634) = lu(634) - lu(313) * lu(615)
         lu(635) = lu(635) - lu(314) * lu(615)
         lu(637) = lu(637) - lu(315) * lu(615)
         lu(638) = lu(638) - lu(316) * lu(615)
         lu(640) = lu(640) - lu(317) * lu(615)
         lu(668) = lu(668) - lu(307) * lu(667)
         lu(671) = lu(671) - lu(308) * lu(667)
         lu(673) = lu(673) - lu(309) * lu(667)
         lu(676) = lu(676) - lu(310) * lu(667)
         lu(677) = lu(677) - lu(311) * lu(667)
         lu(683) = lu(683) - lu(312) * lu(667)
         lu(684) = lu(684) - lu(313) * lu(667)
         lu(685) = lu(685) - lu(314) * lu(667)
         lu(687) = lu(687) - lu(315) * lu(667)
         lu(688) = lu(688) - lu(316) * lu(667)
         lu(690) = lu(690) - lu(317) * lu(667)
         lu(321) = 1._r8 / lu(321)
         lu(322) = lu(322) * lu(321)
         lu(323) = lu(323) * lu(321)
         lu(324) = lu(324) * lu(321)
         lu(325) = lu(325) * lu(321)
         lu(326) = lu(326) * lu(321)
         lu(327) = lu(327) * lu(321)
         lu(328) = lu(328) * lu(321)
         lu(546) = lu(546) - lu(322) * lu(543)
         lu(552) = lu(552) - lu(323) * lu(543)
         lu(558) = lu(558) - lu(324) * lu(543)
         lu(559) = lu(559) - lu(325) * lu(543)
         lu(560) = lu(560) - lu(326) * lu(543)
         lu(562) = lu(562) - lu(327) * lu(543)
         lu(563) = lu(563) - lu(328) * lu(543)
         lu(621) = lu(621) - lu(322) * lu(616)
         lu(627) = lu(627) - lu(323) * lu(616)
         lu(633) = lu(633) - lu(324) * lu(616)
         lu(634) = lu(634) - lu(325) * lu(616)
         lu(635) = lu(635) - lu(326) * lu(616)
         lu(637) = lu(637) - lu(327) * lu(616)
         lu(638) = lu(638) - lu(328) * lu(616)
         lu(671) = lu(671) - lu(322) * lu(668)
         lu(677) = lu(677) - lu(323) * lu(668)
         lu(683) = lu(683) - lu(324) * lu(668)
         lu(684) = lu(684) - lu(325) * lu(668)
         lu(685) = lu(685) - lu(326) * lu(668)
         lu(687) = lu(687) - lu(327) * lu(668)
         lu(688) = lu(688) - lu(328) * lu(668)
         lu(731) = lu(731) - lu(322) * lu(729)
         lu(736) = lu(736) - lu(323) * lu(729)
         lu(742) = lu(742) - lu(324) * lu(729)
         lu(743) = lu(743) - lu(325) * lu(729)
         lu(744) = lu(744) - lu(326) * lu(729)
         lu(746) = lu(746) - lu(327) * lu(729)
         lu(747) = lu(747) - lu(328) * lu(729)
         lu(331) = 1._r8 / lu(331)
         lu(332) = lu(332) * lu(331)
         lu(333) = lu(333) * lu(331)
         lu(334) = lu(334) * lu(331)
         lu(335) = lu(335) * lu(331)
         lu(336) = lu(336) * lu(331)
         lu(337) = lu(337) * lu(331)
         lu(338) = lu(338) * lu(331)
         lu(339) = lu(339) * lu(331)
         lu(340) = lu(340) * lu(331)
         lu(546) = lu(546) - lu(332) * lu(544)
         lu(548) = lu(548) - lu(333) * lu(544)
         lu(552) = lu(552) - lu(334) * lu(544)
         lu(558) = lu(558) - lu(335) * lu(544)
         lu(559) = lu(559) - lu(336) * lu(544)
         lu(560) = lu(560) - lu(337) * lu(544)
         lu(562) = lu(562) - lu(338) * lu(544)
         lu(563) = lu(563) - lu(339) * lu(544)
         lu(565) = lu(565) - lu(340) * lu(544)
         lu(621) = lu(621) - lu(332) * lu(617)
         lu(623) = lu(623) - lu(333) * lu(617)
         lu(627) = lu(627) - lu(334) * lu(617)
         lu(633) = lu(633) - lu(335) * lu(617)
         lu(634) = lu(634) - lu(336) * lu(617)
         lu(635) = lu(635) - lu(337) * lu(617)
         lu(637) = lu(637) - lu(338) * lu(617)
         lu(638) = lu(638) - lu(339) * lu(617)
         lu(640) = lu(640) - lu(340) * lu(617)
         lu(671) = lu(671) - lu(332) * lu(669)
         lu(673) = lu(673) - lu(333) * lu(669)
         lu(677) = lu(677) - lu(334) * lu(669)
         lu(683) = lu(683) - lu(335) * lu(669)
         lu(684) = lu(684) - lu(336) * lu(669)
         lu(685) = lu(685) - lu(337) * lu(669)
         lu(687) = lu(687) - lu(338) * lu(669)
         lu(688) = lu(688) - lu(339) * lu(669)
         lu(690) = lu(690) - lu(340) * lu(669)
         lu(731) = lu(731) - lu(332) * lu(730)
         lu(732) = lu(732) - lu(333) * lu(730)
         lu(736) = lu(736) - lu(334) * lu(730)
         lu(742) = lu(742) - lu(335) * lu(730)
         lu(743) = lu(743) - lu(336) * lu(730)
         lu(744) = lu(744) - lu(337) * lu(730)
         lu(746) = lu(746) - lu(338) * lu(730)
         lu(747) = lu(747) - lu(339) * lu(730)
         lu(749) = lu(749) - lu(340) * lu(730)
      end subroutine lu_fac07
      subroutine lu_fac08( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(342) = 1._r8 / lu(342)
         lu(343) = lu(343) * lu(342)
         lu(344) = lu(344) * lu(342)
         lu(345) = lu(345) * lu(342)
         lu(346) = lu(346) * lu(342)
         lu(347) = lu(347) * lu(342)
         lu(348) = lu(348) * lu(342)
         lu(355) = lu(355) - lu(343) * lu(351)
         lu(356) = - lu(344) * lu(351)
         lu(361) = lu(361) - lu(345) * lu(351)
         lu(362) = lu(362) - lu(346) * lu(351)
         lu(365) = lu(365) - lu(347) * lu(351)
         lu(366) = lu(366) - lu(348) * lu(351)
         lu(480) = lu(480) - lu(343) * lu(477)
         lu(481) = - lu(344) * lu(477)
         lu(485) = lu(485) - lu(345) * lu(477)
         lu(486) = lu(486) - lu(346) * lu(477)
         lu(488) = lu(488) - lu(347) * lu(477)
         lu(489) = lu(489) - lu(348) * lu(477)
         lu(497) = lu(497) - lu(343) * lu(494)
         lu(498) = lu(498) - lu(344) * lu(494)
         lu(504) = lu(504) - lu(345) * lu(494)
         lu(505) = lu(505) - lu(346) * lu(494)
         lu(508) = lu(508) - lu(347) * lu(494)
         lu(509) = lu(509) - lu(348) * lu(494)
         lu(552) = lu(552) - lu(343) * lu(545)
         lu(553) = lu(553) - lu(344) * lu(545)
         lu(559) = lu(559) - lu(345) * lu(545)
         lu(560) = lu(560) - lu(346) * lu(545)
         lu(563) = lu(563) - lu(347) * lu(545)
         lu(564) = lu(564) - lu(348) * lu(545)
         lu(627) = lu(627) - lu(343) * lu(618)
         lu(628) = lu(628) - lu(344) * lu(618)
         lu(634) = lu(634) - lu(345) * lu(618)
         lu(635) = lu(635) - lu(346) * lu(618)
         lu(638) = lu(638) - lu(347) * lu(618)
         lu(639) = lu(639) - lu(348) * lu(618)
         lu(677) = lu(677) - lu(343) * lu(670)
         lu(678) = lu(678) - lu(344) * lu(670)
         lu(684) = lu(684) - lu(345) * lu(670)
         lu(685) = lu(685) - lu(346) * lu(670)
         lu(688) = lu(688) - lu(347) * lu(670)
         lu(689) = lu(689) - lu(348) * lu(670)
         lu(789) = lu(789) - lu(343) * lu(780)
         lu(790) = lu(790) - lu(344) * lu(780)
         lu(796) = lu(796) - lu(345) * lu(780)
         lu(797) = lu(797) - lu(346) * lu(780)
         lu(800) = lu(800) - lu(347) * lu(780)
         lu(801) = lu(801) - lu(348) * lu(780)
         lu(352) = 1._r8 / lu(352)
         lu(353) = lu(353) * lu(352)
         lu(354) = lu(354) * lu(352)
         lu(355) = lu(355) * lu(352)
         lu(356) = lu(356) * lu(352)
         lu(357) = lu(357) * lu(352)
         lu(358) = lu(358) * lu(352)
         lu(359) = lu(359) * lu(352)
         lu(360) = lu(360) * lu(352)
         lu(361) = lu(361) * lu(352)
         lu(362) = lu(362) * lu(352)
         lu(363) = lu(363) * lu(352)
         lu(364) = lu(364) * lu(352)
         lu(365) = lu(365) * lu(352)
         lu(366) = lu(366) * lu(352)
         lu(620) = lu(620) - lu(353) * lu(619)
         lu(625) = lu(625) - lu(354) * lu(619)
         lu(627) = lu(627) - lu(355) * lu(619)
         lu(628) = lu(628) - lu(356) * lu(619)
         lu(629) = lu(629) - lu(357) * lu(619)
         lu(631) = lu(631) - lu(358) * lu(619)
         lu(632) = lu(632) - lu(359) * lu(619)
         lu(633) = lu(633) - lu(360) * lu(619)
         lu(634) = lu(634) - lu(361) * lu(619)
         lu(635) = lu(635) - lu(362) * lu(619)
         lu(636) = lu(636) - lu(363) * lu(619)
         lu(637) = lu(637) - lu(364) * lu(619)
         lu(638) = lu(638) - lu(365) * lu(619)
         lu(639) = lu(639) - lu(366) * lu(619)
         lu(700) = lu(700) - lu(353) * lu(699)
         lu(705) = lu(705) - lu(354) * lu(699)
         lu(707) = lu(707) - lu(355) * lu(699)
         lu(708) = - lu(356) * lu(699)
         lu(709) = lu(709) - lu(357) * lu(699)
         lu(711) = - lu(358) * lu(699)
         lu(712) = lu(712) - lu(359) * lu(699)
         lu(713) = lu(713) - lu(360) * lu(699)
         lu(714) = lu(714) - lu(361) * lu(699)
         lu(715) = lu(715) - lu(362) * lu(699)
         lu(716) = lu(716) - lu(363) * lu(699)
         lu(717) = lu(717) - lu(364) * lu(699)
         lu(718) = lu(718) - lu(365) * lu(699)
         lu(719) = lu(719) - lu(366) * lu(699)
         lu(782) = lu(782) - lu(353) * lu(781)
         lu(787) = lu(787) - lu(354) * lu(781)
         lu(789) = lu(789) - lu(355) * lu(781)
         lu(790) = lu(790) - lu(356) * lu(781)
         lu(791) = lu(791) - lu(357) * lu(781)
         lu(793) = lu(793) - lu(358) * lu(781)
         lu(794) = lu(794) - lu(359) * lu(781)
         lu(795) = lu(795) - lu(360) * lu(781)
         lu(796) = lu(796) - lu(361) * lu(781)
         lu(797) = lu(797) - lu(362) * lu(781)
         lu(798) = lu(798) - lu(363) * lu(781)
         lu(799) = - lu(364) * lu(781)
         lu(800) = lu(800) - lu(365) * lu(781)
         lu(801) = lu(801) - lu(366) * lu(781)
         lu(370) = 1._r8 / lu(370)
         lu(371) = lu(371) * lu(370)
         lu(372) = lu(372) * lu(370)
         lu(373) = lu(373) * lu(370)
         lu(374) = lu(374) * lu(370)
         lu(375) = lu(375) * lu(370)
         lu(376) = lu(376) * lu(370)
         lu(377) = lu(377) * lu(370)
         lu(378) = lu(378) * lu(370)
         lu(379) = lu(379) * lu(370)
         lu(380) = lu(380) * lu(370)
         lu(381) = lu(381) * lu(370)
         lu(460) = lu(460) - lu(371) * lu(459)
         lu(461) = - lu(372) * lu(459)
         lu(462) = lu(462) - lu(373) * lu(459)
         lu(463) = lu(463) - lu(374) * lu(459)
         lu(466) = - lu(375) * lu(459)
         lu(467) = lu(467) - lu(376) * lu(459)
         lu(468) = lu(468) - lu(377) * lu(459)
         lu(469) = lu(469) - lu(378) * lu(459)
         lu(470) = lu(470) - lu(379) * lu(459)
         lu(471) = - lu(380) * lu(459)
         lu(472) = - lu(381) * lu(459)
         lu(621) = lu(621) - lu(371) * lu(620)
         lu(624) = lu(624) - lu(372) * lu(620)
         lu(626) = lu(626) - lu(373) * lu(620)
         lu(627) = lu(627) - lu(374) * lu(620)
         lu(633) = lu(633) - lu(375) * lu(620)
         lu(634) = lu(634) - lu(376) * lu(620)
         lu(635) = lu(635) - lu(377) * lu(620)
         lu(636) = lu(636) - lu(378) * lu(620)
         lu(637) = lu(637) - lu(379) * lu(620)
         lu(638) = lu(638) - lu(380) * lu(620)
         lu(639) = lu(639) - lu(381) * lu(620)
         lu(701) = lu(701) - lu(371) * lu(700)
         lu(704) = - lu(372) * lu(700)
         lu(706) = lu(706) - lu(373) * lu(700)
         lu(707) = lu(707) - lu(374) * lu(700)
         lu(713) = lu(713) - lu(375) * lu(700)
         lu(714) = lu(714) - lu(376) * lu(700)
         lu(715) = lu(715) - lu(377) * lu(700)
         lu(716) = lu(716) - lu(378) * lu(700)
         lu(717) = lu(717) - lu(379) * lu(700)
         lu(718) = lu(718) - lu(380) * lu(700)
         lu(719) = lu(719) - lu(381) * lu(700)
         lu(783) = lu(783) - lu(371) * lu(782)
         lu(786) = lu(786) - lu(372) * lu(782)
         lu(788) = lu(788) - lu(373) * lu(782)
         lu(789) = lu(789) - lu(374) * lu(782)
         lu(795) = lu(795) - lu(375) * lu(782)
         lu(796) = lu(796) - lu(376) * lu(782)
         lu(797) = lu(797) - lu(377) * lu(782)
         lu(798) = lu(798) - lu(378) * lu(782)
         lu(799) = lu(799) - lu(379) * lu(782)
         lu(800) = lu(800) - lu(380) * lu(782)
         lu(801) = lu(801) - lu(381) * lu(782)
         lu(383) = 1._r8 / lu(383)
         lu(384) = lu(384) * lu(383)
         lu(385) = lu(385) * lu(383)
         lu(386) = lu(386) * lu(383)
         lu(387) = lu(387) * lu(383)
         lu(388) = lu(388) * lu(383)
         lu(389) = lu(389) * lu(383)
         lu(467) = lu(467) - lu(384) * lu(460)
         lu(468) = lu(468) - lu(385) * lu(460)
         lu(470) = lu(470) - lu(386) * lu(460)
         lu(471) = lu(471) - lu(387) * lu(460)
         lu(472) = lu(472) - lu(388) * lu(460)
         lu(473) = lu(473) - lu(389) * lu(460)
         lu(559) = lu(559) - lu(384) * lu(546)
         lu(560) = lu(560) - lu(385) * lu(546)
         lu(562) = lu(562) - lu(386) * lu(546)
         lu(563) = lu(563) - lu(387) * lu(546)
         lu(564) = lu(564) - lu(388) * lu(546)
         lu(565) = lu(565) - lu(389) * lu(546)
         lu(634) = lu(634) - lu(384) * lu(621)
         lu(635) = lu(635) - lu(385) * lu(621)
         lu(637) = lu(637) - lu(386) * lu(621)
         lu(638) = lu(638) - lu(387) * lu(621)
         lu(639) = lu(639) - lu(388) * lu(621)
         lu(640) = lu(640) - lu(389) * lu(621)
         lu(684) = lu(684) - lu(384) * lu(671)
         lu(685) = lu(685) - lu(385) * lu(671)
         lu(687) = lu(687) - lu(386) * lu(671)
         lu(688) = lu(688) - lu(387) * lu(671)
         lu(689) = lu(689) - lu(388) * lu(671)
         lu(690) = lu(690) - lu(389) * lu(671)
         lu(714) = lu(714) - lu(384) * lu(701)
         lu(715) = lu(715) - lu(385) * lu(701)
         lu(717) = lu(717) - lu(386) * lu(701)
         lu(718) = lu(718) - lu(387) * lu(701)
         lu(719) = lu(719) - lu(388) * lu(701)
         lu(720) = - lu(389) * lu(701)
         lu(743) = lu(743) - lu(384) * lu(731)
         lu(744) = lu(744) - lu(385) * lu(731)
         lu(746) = lu(746) - lu(386) * lu(731)
         lu(747) = lu(747) - lu(387) * lu(731)
         lu(748) = - lu(388) * lu(731)
         lu(749) = lu(749) - lu(389) * lu(731)
         lu(796) = lu(796) - lu(384) * lu(783)
         lu(797) = lu(797) - lu(385) * lu(783)
         lu(799) = lu(799) - lu(386) * lu(783)
         lu(800) = lu(800) - lu(387) * lu(783)
         lu(801) = lu(801) - lu(388) * lu(783)
         lu(802) = lu(802) - lu(389) * lu(783)
      end subroutine lu_fac08
      subroutine lu_fac09( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(392) = 1._r8 / lu(392)
         lu(393) = lu(393) * lu(392)
         lu(394) = lu(394) * lu(392)
         lu(395) = lu(395) * lu(392)
         lu(396) = lu(396) * lu(392)
         lu(397) = lu(397) * lu(392)
         lu(398) = lu(398) * lu(392)
         lu(399) = lu(399) * lu(392)
         lu(400) = lu(400) * lu(392)
         lu(401) = lu(401) * lu(392)
         lu(548) = lu(548) - lu(393) * lu(547)
         lu(550) = lu(550) - lu(394) * lu(547)
         lu(554) = lu(554) - lu(395) * lu(547)
         lu(558) = lu(558) - lu(396) * lu(547)
         lu(559) = lu(559) - lu(397) * lu(547)
         lu(560) = lu(560) - lu(398) * lu(547)
         lu(562) = lu(562) - lu(399) * lu(547)
         lu(563) = lu(563) - lu(400) * lu(547)
         lu(565) = lu(565) - lu(401) * lu(547)
         lu(623) = lu(623) - lu(393) * lu(622)
         lu(625) = lu(625) - lu(394) * lu(622)
         lu(629) = lu(629) - lu(395) * lu(622)
         lu(633) = lu(633) - lu(396) * lu(622)
         lu(634) = lu(634) - lu(397) * lu(622)
         lu(635) = lu(635) - lu(398) * lu(622)
         lu(637) = lu(637) - lu(399) * lu(622)
         lu(638) = lu(638) - lu(400) * lu(622)
         lu(640) = lu(640) - lu(401) * lu(622)
         lu(673) = lu(673) - lu(393) * lu(672)
         lu(675) = lu(675) - lu(394) * lu(672)
         lu(679) = lu(679) - lu(395) * lu(672)
         lu(683) = lu(683) - lu(396) * lu(672)
         lu(684) = lu(684) - lu(397) * lu(672)
         lu(685) = lu(685) - lu(398) * lu(672)
         lu(687) = lu(687) - lu(399) * lu(672)
         lu(688) = lu(688) - lu(400) * lu(672)
         lu(690) = lu(690) - lu(401) * lu(672)
         lu(703) = - lu(393) * lu(702)
         lu(705) = lu(705) - lu(394) * lu(702)
         lu(709) = lu(709) - lu(395) * lu(702)
         lu(713) = lu(713) - lu(396) * lu(702)
         lu(714) = lu(714) - lu(397) * lu(702)
         lu(715) = lu(715) - lu(398) * lu(702)
         lu(717) = lu(717) - lu(399) * lu(702)
         lu(718) = lu(718) - lu(400) * lu(702)
         lu(720) = lu(720) - lu(401) * lu(702)
         lu(785) = - lu(393) * lu(784)
         lu(787) = lu(787) - lu(394) * lu(784)
         lu(791) = lu(791) - lu(395) * lu(784)
         lu(795) = lu(795) - lu(396) * lu(784)
         lu(796) = lu(796) - lu(397) * lu(784)
         lu(797) = lu(797) - lu(398) * lu(784)
         lu(799) = lu(799) - lu(399) * lu(784)
         lu(800) = lu(800) - lu(400) * lu(784)
         lu(802) = lu(802) - lu(401) * lu(784)
         lu(404) = 1._r8 / lu(404)
         lu(405) = lu(405) * lu(404)
         lu(406) = lu(406) * lu(404)
         lu(407) = lu(407) * lu(404)
         lu(408) = lu(408) * lu(404)
         lu(409) = lu(409) * lu(404)
         lu(410) = lu(410) * lu(404)
         lu(411) = lu(411) * lu(404)
         lu(412) = lu(412) * lu(404)
         lu(413) = lu(413) * lu(404)
         lu(549) = lu(549) - lu(405) * lu(548)
         lu(551) = lu(551) - lu(406) * lu(548)
         lu(552) = lu(552) - lu(407) * lu(548)
         lu(558) = lu(558) - lu(408) * lu(548)
         lu(559) = lu(559) - lu(409) * lu(548)
         lu(560) = lu(560) - lu(410) * lu(548)
         lu(562) = lu(562) - lu(411) * lu(548)
         lu(563) = lu(563) - lu(412) * lu(548)
         lu(565) = lu(565) - lu(413) * lu(548)
         lu(624) = lu(624) - lu(405) * lu(623)
         lu(626) = lu(626) - lu(406) * lu(623)
         lu(627) = lu(627) - lu(407) * lu(623)
         lu(633) = lu(633) - lu(408) * lu(623)
         lu(634) = lu(634) - lu(409) * lu(623)
         lu(635) = lu(635) - lu(410) * lu(623)
         lu(637) = lu(637) - lu(411) * lu(623)
         lu(638) = lu(638) - lu(412) * lu(623)
         lu(640) = lu(640) - lu(413) * lu(623)
         lu(674) = lu(674) - lu(405) * lu(673)
         lu(676) = lu(676) - lu(406) * lu(673)
         lu(677) = lu(677) - lu(407) * lu(673)
         lu(683) = lu(683) - lu(408) * lu(673)
         lu(684) = lu(684) - lu(409) * lu(673)
         lu(685) = lu(685) - lu(410) * lu(673)
         lu(687) = lu(687) - lu(411) * lu(673)
         lu(688) = lu(688) - lu(412) * lu(673)
         lu(690) = lu(690) - lu(413) * lu(673)
         lu(704) = lu(704) - lu(405) * lu(703)
         lu(706) = lu(706) - lu(406) * lu(703)
         lu(707) = lu(707) - lu(407) * lu(703)
         lu(713) = lu(713) - lu(408) * lu(703)
         lu(714) = lu(714) - lu(409) * lu(703)
         lu(715) = lu(715) - lu(410) * lu(703)
         lu(717) = lu(717) - lu(411) * lu(703)
         lu(718) = lu(718) - lu(412) * lu(703)
         lu(720) = lu(720) - lu(413) * lu(703)
         lu(733) = lu(733) - lu(405) * lu(732)
         lu(735) = lu(735) - lu(406) * lu(732)
         lu(736) = lu(736) - lu(407) * lu(732)
         lu(742) = lu(742) - lu(408) * lu(732)
         lu(743) = lu(743) - lu(409) * lu(732)
         lu(744) = lu(744) - lu(410) * lu(732)
         lu(746) = lu(746) - lu(411) * lu(732)
         lu(747) = lu(747) - lu(412) * lu(732)
         lu(749) = lu(749) - lu(413) * lu(732)
         lu(786) = lu(786) - lu(405) * lu(785)
         lu(788) = lu(788) - lu(406) * lu(785)
         lu(789) = lu(789) - lu(407) * lu(785)
         lu(795) = lu(795) - lu(408) * lu(785)
         lu(796) = lu(796) - lu(409) * lu(785)
         lu(797) = lu(797) - lu(410) * lu(785)
         lu(799) = lu(799) - lu(411) * lu(785)
         lu(800) = lu(800) - lu(412) * lu(785)
         lu(802) = lu(802) - lu(413) * lu(785)
         lu(414) = 1._r8 / lu(414)
         lu(415) = lu(415) * lu(414)
         lu(416) = lu(416) * lu(414)
         lu(417) = lu(417) * lu(414)
         lu(418) = lu(418) * lu(414)
         lu(419) = lu(419) * lu(414)
         lu(449) = lu(449) - lu(415) * lu(448)
         lu(450) = lu(450) - lu(416) * lu(448)
         lu(453) = lu(453) - lu(417) * lu(448)
         lu(454) = lu(454) - lu(418) * lu(448)
         lu(458) = lu(458) - lu(419) * lu(448)
         lu(462) = lu(462) - lu(415) * lu(461)
         lu(463) = lu(463) - lu(416) * lu(461)
         lu(467) = lu(467) - lu(417) * lu(461)
         lu(468) = lu(468) - lu(418) * lu(461)
         lu(473) = lu(473) - lu(419) * lu(461)
         lu(479) = lu(479) - lu(415) * lu(478)
         lu(480) = lu(480) - lu(416) * lu(478)
         lu(485) = lu(485) - lu(417) * lu(478)
         lu(486) = lu(486) - lu(418) * lu(478)
         lu(490) = lu(490) - lu(419) * lu(478)
         lu(515) = - lu(415) * lu(514)
         lu(516) = lu(516) - lu(416) * lu(514)
         lu(519) = lu(519) - lu(417) * lu(514)
         lu(520) = lu(520) - lu(418) * lu(514)
         lu(525) = lu(525) - lu(419) * lu(514)
         lu(551) = lu(551) - lu(415) * lu(549)
         lu(552) = lu(552) - lu(416) * lu(549)
         lu(559) = lu(559) - lu(417) * lu(549)
         lu(560) = lu(560) - lu(418) * lu(549)
         lu(565) = lu(565) - lu(419) * lu(549)
         lu(626) = lu(626) - lu(415) * lu(624)
         lu(627) = lu(627) - lu(416) * lu(624)
         lu(634) = lu(634) - lu(417) * lu(624)
         lu(635) = lu(635) - lu(418) * lu(624)
         lu(640) = lu(640) - lu(419) * lu(624)
         lu(676) = lu(676) - lu(415) * lu(674)
         lu(677) = lu(677) - lu(416) * lu(674)
         lu(684) = lu(684) - lu(417) * lu(674)
         lu(685) = lu(685) - lu(418) * lu(674)
         lu(690) = lu(690) - lu(419) * lu(674)
         lu(706) = lu(706) - lu(415) * lu(704)
         lu(707) = lu(707) - lu(416) * lu(704)
         lu(714) = lu(714) - lu(417) * lu(704)
         lu(715) = lu(715) - lu(418) * lu(704)
         lu(720) = lu(720) - lu(419) * lu(704)
         lu(735) = lu(735) - lu(415) * lu(733)
         lu(736) = lu(736) - lu(416) * lu(733)
         lu(743) = lu(743) - lu(417) * lu(733)
         lu(744) = lu(744) - lu(418) * lu(733)
         lu(749) = lu(749) - lu(419) * lu(733)
         lu(759) = lu(759) - lu(415) * lu(758)
         lu(760) = lu(760) - lu(416) * lu(758)
         lu(763) = lu(763) - lu(417) * lu(758)
         lu(764) = lu(764) - lu(418) * lu(758)
         lu(769) = lu(769) - lu(419) * lu(758)
         lu(788) = lu(788) - lu(415) * lu(786)
         lu(789) = lu(789) - lu(416) * lu(786)
         lu(796) = lu(796) - lu(417) * lu(786)
         lu(797) = lu(797) - lu(418) * lu(786)
         lu(802) = lu(802) - lu(419) * lu(786)
         lu(810) = lu(810) - lu(415) * lu(808)
         lu(811) = lu(811) - lu(416) * lu(808)
         lu(818) = lu(818) - lu(417) * lu(808)
         lu(819) = lu(819) - lu(418) * lu(808)
         lu(824) = lu(824) - lu(419) * lu(808)
         lu(420) = 1._r8 / lu(420)
         lu(421) = lu(421) * lu(420)
         lu(422) = lu(422) * lu(420)
         lu(423) = lu(423) * lu(420)
         lu(424) = lu(424) * lu(420)
         lu(425) = lu(425) * lu(420)
         lu(426) = lu(426) * lu(420)
         lu(427) = lu(427) * lu(420)
         lu(428) = lu(428) * lu(420)
         lu(496) = - lu(421) * lu(495)
         lu(497) = lu(497) - lu(422) * lu(495)
         lu(500) = - lu(423) * lu(495)
         lu(502) = - lu(424) * lu(495)
         lu(504) = lu(504) - lu(425) * lu(495)
         lu(505) = lu(505) - lu(426) * lu(495)
         lu(506) = - lu(427) * lu(495)
         lu(510) = lu(510) - lu(428) * lu(495)
         lu(551) = lu(551) - lu(421) * lu(550)
         lu(552) = lu(552) - lu(422) * lu(550)
         lu(555) = lu(555) - lu(423) * lu(550)
         lu(557) = lu(557) - lu(424) * lu(550)
         lu(559) = lu(559) - lu(425) * lu(550)
         lu(560) = lu(560) - lu(426) * lu(550)
         lu(561) = lu(561) - lu(427) * lu(550)
         lu(565) = lu(565) - lu(428) * lu(550)
         lu(626) = lu(626) - lu(421) * lu(625)
         lu(627) = lu(627) - lu(422) * lu(625)
         lu(630) = lu(630) - lu(423) * lu(625)
         lu(632) = lu(632) - lu(424) * lu(625)
         lu(634) = lu(634) - lu(425) * lu(625)
         lu(635) = lu(635) - lu(426) * lu(625)
         lu(636) = lu(636) - lu(427) * lu(625)
         lu(640) = lu(640) - lu(428) * lu(625)
         lu(676) = lu(676) - lu(421) * lu(675)
         lu(677) = lu(677) - lu(422) * lu(675)
         lu(680) = lu(680) - lu(423) * lu(675)
         lu(682) = lu(682) - lu(424) * lu(675)
         lu(684) = lu(684) - lu(425) * lu(675)
         lu(685) = lu(685) - lu(426) * lu(675)
         lu(686) = lu(686) - lu(427) * lu(675)
         lu(690) = lu(690) - lu(428) * lu(675)
         lu(706) = lu(706) - lu(421) * lu(705)
         lu(707) = lu(707) - lu(422) * lu(705)
         lu(710) = - lu(423) * lu(705)
         lu(712) = lu(712) - lu(424) * lu(705)
         lu(714) = lu(714) - lu(425) * lu(705)
         lu(715) = lu(715) - lu(426) * lu(705)
         lu(716) = lu(716) - lu(427) * lu(705)
         lu(720) = lu(720) - lu(428) * lu(705)
         lu(735) = lu(735) - lu(421) * lu(734)
         lu(736) = lu(736) - lu(422) * lu(734)
         lu(739) = lu(739) - lu(423) * lu(734)
         lu(741) = lu(741) - lu(424) * lu(734)
         lu(743) = lu(743) - lu(425) * lu(734)
         lu(744) = lu(744) - lu(426) * lu(734)
         lu(745) = - lu(427) * lu(734)
         lu(749) = lu(749) - lu(428) * lu(734)
         lu(788) = lu(788) - lu(421) * lu(787)
         lu(789) = lu(789) - lu(422) * lu(787)
         lu(792) = lu(792) - lu(423) * lu(787)
         lu(794) = lu(794) - lu(424) * lu(787)
         lu(796) = lu(796) - lu(425) * lu(787)
         lu(797) = lu(797) - lu(426) * lu(787)
         lu(798) = lu(798) - lu(427) * lu(787)
         lu(802) = lu(802) - lu(428) * lu(787)
         lu(810) = lu(810) - lu(421) * lu(809)
         lu(811) = lu(811) - lu(422) * lu(809)
         lu(814) = lu(814) - lu(423) * lu(809)
         lu(816) = lu(816) - lu(424) * lu(809)
         lu(818) = lu(818) - lu(425) * lu(809)
         lu(819) = lu(819) - lu(426) * lu(809)
         lu(820) = lu(820) - lu(427) * lu(809)
         lu(824) = lu(824) - lu(428) * lu(809)
         lu(430) = 1._r8 / lu(430)
         lu(431) = lu(431) * lu(430)
         lu(432) = lu(432) * lu(430)
         lu(433) = lu(433) * lu(430)
         lu(434) = lu(434) * lu(430)
         lu(435) = lu(435) * lu(430)
         lu(453) = lu(453) - lu(431) * lu(449)
         lu(454) = lu(454) - lu(432) * lu(449)
         lu(456) = lu(456) - lu(433) * lu(449)
         lu(457) = lu(457) - lu(434) * lu(449)
         lu(458) = lu(458) - lu(435) * lu(449)
         lu(467) = lu(467) - lu(431) * lu(462)
         lu(468) = lu(468) - lu(432) * lu(462)
         lu(471) = lu(471) - lu(433) * lu(462)
         lu(472) = lu(472) - lu(434) * lu(462)
         lu(473) = lu(473) - lu(435) * lu(462)
         lu(485) = lu(485) - lu(431) * lu(479)
         lu(486) = lu(486) - lu(432) * lu(479)
         lu(488) = lu(488) - lu(433) * lu(479)
         lu(489) = lu(489) - lu(434) * lu(479)
         lu(490) = lu(490) - lu(435) * lu(479)
         lu(504) = lu(504) - lu(431) * lu(496)
         lu(505) = lu(505) - lu(432) * lu(496)
         lu(508) = lu(508) - lu(433) * lu(496)
         lu(509) = lu(509) - lu(434) * lu(496)
         lu(510) = lu(510) - lu(435) * lu(496)
         lu(519) = lu(519) - lu(431) * lu(515)
         lu(520) = lu(520) - lu(432) * lu(515)
         lu(523) = lu(523) - lu(433) * lu(515)
         lu(524) = lu(524) - lu(434) * lu(515)
         lu(525) = lu(525) - lu(435) * lu(515)
         lu(559) = lu(559) - lu(431) * lu(551)
         lu(560) = lu(560) - lu(432) * lu(551)
         lu(563) = lu(563) - lu(433) * lu(551)
         lu(564) = lu(564) - lu(434) * lu(551)
         lu(565) = lu(565) - lu(435) * lu(551)
         lu(634) = lu(634) - lu(431) * lu(626)
         lu(635) = lu(635) - lu(432) * lu(626)
         lu(638) = lu(638) - lu(433) * lu(626)
         lu(639) = lu(639) - lu(434) * lu(626)
         lu(640) = lu(640) - lu(435) * lu(626)
         lu(684) = lu(684) - lu(431) * lu(676)
         lu(685) = lu(685) - lu(432) * lu(676)
         lu(688) = lu(688) - lu(433) * lu(676)
         lu(689) = lu(689) - lu(434) * lu(676)
         lu(690) = lu(690) - lu(435) * lu(676)
         lu(714) = lu(714) - lu(431) * lu(706)
         lu(715) = lu(715) - lu(432) * lu(706)
         lu(718) = lu(718) - lu(433) * lu(706)
         lu(719) = lu(719) - lu(434) * lu(706)
         lu(720) = lu(720) - lu(435) * lu(706)
         lu(743) = lu(743) - lu(431) * lu(735)
         lu(744) = lu(744) - lu(432) * lu(735)
         lu(747) = lu(747) - lu(433) * lu(735)
         lu(748) = lu(748) - lu(434) * lu(735)
         lu(749) = lu(749) - lu(435) * lu(735)
         lu(763) = lu(763) - lu(431) * lu(759)
         lu(764) = lu(764) - lu(432) * lu(759)
         lu(767) = lu(767) - lu(433) * lu(759)
         lu(768) = lu(768) - lu(434) * lu(759)
         lu(769) = lu(769) - lu(435) * lu(759)
         lu(796) = lu(796) - lu(431) * lu(788)
         lu(797) = lu(797) - lu(432) * lu(788)
         lu(800) = lu(800) - lu(433) * lu(788)
         lu(801) = lu(801) - lu(434) * lu(788)
         lu(802) = lu(802) - lu(435) * lu(788)
         lu(818) = lu(818) - lu(431) * lu(810)
         lu(819) = lu(819) - lu(432) * lu(810)
         lu(822) = lu(822) - lu(433) * lu(810)
         lu(823) = lu(823) - lu(434) * lu(810)
         lu(824) = lu(824) - lu(435) * lu(810)
         lu(438) = 1._r8 / lu(438)
         lu(439) = lu(439) * lu(438)
         lu(440) = lu(440) * lu(438)
         lu(441) = lu(441) * lu(438)
         lu(442) = lu(442) * lu(438)
         lu(443) = lu(443) * lu(438)
         lu(452) = lu(452) - lu(439) * lu(450)
         lu(453) = lu(453) - lu(440) * lu(450)
         lu(454) = lu(454) - lu(441) * lu(450)
         lu(456) = lu(456) - lu(442) * lu(450)
         lu(457) = lu(457) - lu(443) * lu(450)
         lu(466) = lu(466) - lu(439) * lu(463)
         lu(467) = lu(467) - lu(440) * lu(463)
         lu(468) = lu(468) - lu(441) * lu(463)
         lu(471) = lu(471) - lu(442) * lu(463)
         lu(472) = lu(472) - lu(443) * lu(463)
         lu(484) = lu(484) - lu(439) * lu(480)
         lu(485) = lu(485) - lu(440) * lu(480)
         lu(486) = lu(486) - lu(441) * lu(480)
         lu(488) = lu(488) - lu(442) * lu(480)
         lu(489) = lu(489) - lu(443) * lu(480)
         lu(503) = lu(503) - lu(439) * lu(497)
         lu(504) = lu(504) - lu(440) * lu(497)
         lu(505) = lu(505) - lu(441) * lu(497)
         lu(508) = lu(508) - lu(442) * lu(497)
         lu(509) = lu(509) - lu(443) * lu(497)
         lu(518) = lu(518) - lu(439) * lu(516)
         lu(519) = lu(519) - lu(440) * lu(516)
         lu(520) = lu(520) - lu(441) * lu(516)
         lu(523) = lu(523) - lu(442) * lu(516)
         lu(524) = lu(524) - lu(443) * lu(516)
         lu(558) = lu(558) - lu(439) * lu(552)
         lu(559) = lu(559) - lu(440) * lu(552)
         lu(560) = lu(560) - lu(441) * lu(552)
         lu(563) = lu(563) - lu(442) * lu(552)
         lu(564) = lu(564) - lu(443) * lu(552)
         lu(633) = lu(633) - lu(439) * lu(627)
         lu(634) = lu(634) - lu(440) * lu(627)
         lu(635) = lu(635) - lu(441) * lu(627)
         lu(638) = lu(638) - lu(442) * lu(627)
         lu(639) = lu(639) - lu(443) * lu(627)
         lu(683) = lu(683) - lu(439) * lu(677)
         lu(684) = lu(684) - lu(440) * lu(677)
         lu(685) = lu(685) - lu(441) * lu(677)
         lu(688) = lu(688) - lu(442) * lu(677)
         lu(689) = lu(689) - lu(443) * lu(677)
         lu(713) = lu(713) - lu(439) * lu(707)
         lu(714) = lu(714) - lu(440) * lu(707)
         lu(715) = lu(715) - lu(441) * lu(707)
         lu(718) = lu(718) - lu(442) * lu(707)
         lu(719) = lu(719) - lu(443) * lu(707)
         lu(742) = lu(742) - lu(439) * lu(736)
         lu(743) = lu(743) - lu(440) * lu(736)
         lu(744) = lu(744) - lu(441) * lu(736)
         lu(747) = lu(747) - lu(442) * lu(736)
         lu(748) = lu(748) - lu(443) * lu(736)
         lu(762) = lu(762) - lu(439) * lu(760)
         lu(763) = lu(763) - lu(440) * lu(760)
         lu(764) = lu(764) - lu(441) * lu(760)
         lu(767) = lu(767) - lu(442) * lu(760)
         lu(768) = lu(768) - lu(443) * lu(760)
         lu(795) = lu(795) - lu(439) * lu(789)
         lu(796) = lu(796) - lu(440) * lu(789)
         lu(797) = lu(797) - lu(441) * lu(789)
         lu(800) = lu(800) - lu(442) * lu(789)
         lu(801) = lu(801) - lu(443) * lu(789)
         lu(817) = lu(817) - lu(439) * lu(811)
         lu(818) = lu(818) - lu(440) * lu(811)
         lu(819) = lu(819) - lu(441) * lu(811)
         lu(822) = lu(822) - lu(442) * lu(811)
         lu(823) = lu(823) - lu(443) * lu(811)
      end subroutine lu_fac09
      subroutine lu_fac10( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(451) = 1._r8 / lu(451)
         lu(452) = lu(452) * lu(451)
         lu(453) = lu(453) * lu(451)
         lu(454) = lu(454) * lu(451)
         lu(455) = lu(455) * lu(451)
         lu(456) = lu(456) * lu(451)
         lu(457) = lu(457) * lu(451)
         lu(458) = lu(458) * lu(451)
         lu(484) = lu(484) - lu(452) * lu(481)
         lu(485) = lu(485) - lu(453) * lu(481)
         lu(486) = lu(486) - lu(454) * lu(481)
         lu(487) = lu(487) - lu(455) * lu(481)
         lu(488) = lu(488) - lu(456) * lu(481)
         lu(489) = lu(489) - lu(457) * lu(481)
         lu(490) = lu(490) - lu(458) * lu(481)
         lu(503) = lu(503) - lu(452) * lu(498)
         lu(504) = lu(504) - lu(453) * lu(498)
         lu(505) = lu(505) - lu(454) * lu(498)
         lu(507) = lu(507) - lu(455) * lu(498)
         lu(508) = lu(508) - lu(456) * lu(498)
         lu(509) = lu(509) - lu(457) * lu(498)
         lu(510) = lu(510) - lu(458) * lu(498)
         lu(558) = lu(558) - lu(452) * lu(553)
         lu(559) = lu(559) - lu(453) * lu(553)
         lu(560) = lu(560) - lu(454) * lu(553)
         lu(562) = lu(562) - lu(455) * lu(553)
         lu(563) = lu(563) - lu(456) * lu(553)
         lu(564) = lu(564) - lu(457) * lu(553)
         lu(565) = lu(565) - lu(458) * lu(553)
         lu(633) = lu(633) - lu(452) * lu(628)
         lu(634) = lu(634) - lu(453) * lu(628)
         lu(635) = lu(635) - lu(454) * lu(628)
         lu(637) = lu(637) - lu(455) * lu(628)
         lu(638) = lu(638) - lu(456) * lu(628)
         lu(639) = lu(639) - lu(457) * lu(628)
         lu(640) = lu(640) - lu(458) * lu(628)
         lu(683) = lu(683) - lu(452) * lu(678)
         lu(684) = lu(684) - lu(453) * lu(678)
         lu(685) = lu(685) - lu(454) * lu(678)
         lu(687) = lu(687) - lu(455) * lu(678)
         lu(688) = lu(688) - lu(456) * lu(678)
         lu(689) = lu(689) - lu(457) * lu(678)
         lu(690) = lu(690) - lu(458) * lu(678)
         lu(713) = lu(713) - lu(452) * lu(708)
         lu(714) = lu(714) - lu(453) * lu(708)
         lu(715) = lu(715) - lu(454) * lu(708)
         lu(717) = lu(717) - lu(455) * lu(708)
         lu(718) = lu(718) - lu(456) * lu(708)
         lu(719) = lu(719) - lu(457) * lu(708)
         lu(720) = lu(720) - lu(458) * lu(708)
         lu(742) = lu(742) - lu(452) * lu(737)
         lu(743) = lu(743) - lu(453) * lu(737)
         lu(744) = lu(744) - lu(454) * lu(737)
         lu(746) = lu(746) - lu(455) * lu(737)
         lu(747) = lu(747) - lu(456) * lu(737)
         lu(748) = lu(748) - lu(457) * lu(737)
         lu(749) = lu(749) - lu(458) * lu(737)
         lu(795) = lu(795) - lu(452) * lu(790)
         lu(796) = lu(796) - lu(453) * lu(790)
         lu(797) = lu(797) - lu(454) * lu(790)
         lu(799) = lu(799) - lu(455) * lu(790)
         lu(800) = lu(800) - lu(456) * lu(790)
         lu(801) = lu(801) - lu(457) * lu(790)
         lu(802) = lu(802) - lu(458) * lu(790)
         lu(817) = lu(817) - lu(452) * lu(812)
         lu(818) = lu(818) - lu(453) * lu(812)
         lu(819) = lu(819) - lu(454) * lu(812)
         lu(821) = lu(821) - lu(455) * lu(812)
         lu(822) = lu(822) - lu(456) * lu(812)
         lu(823) = lu(823) - lu(457) * lu(812)
         lu(824) = lu(824) - lu(458) * lu(812)
         lu(464) = 1._r8 / lu(464)
         lu(465) = lu(465) * lu(464)
         lu(466) = lu(466) * lu(464)
         lu(467) = lu(467) * lu(464)
         lu(468) = lu(468) * lu(464)
         lu(469) = lu(469) * lu(464)
         lu(470) = lu(470) * lu(464)
         lu(471) = lu(471) * lu(464)
         lu(472) = lu(472) * lu(464)
         lu(473) = lu(473) * lu(464)
         lu(500) = lu(500) - lu(465) * lu(499)
         lu(503) = lu(503) - lu(466) * lu(499)
         lu(504) = lu(504) - lu(467) * lu(499)
         lu(505) = lu(505) - lu(468) * lu(499)
         lu(506) = lu(506) - lu(469) * lu(499)
         lu(507) = lu(507) - lu(470) * lu(499)
         lu(508) = lu(508) - lu(471) * lu(499)
         lu(509) = lu(509) - lu(472) * lu(499)
         lu(510) = lu(510) - lu(473) * lu(499)
         lu(555) = lu(555) - lu(465) * lu(554)
         lu(558) = lu(558) - lu(466) * lu(554)
         lu(559) = lu(559) - lu(467) * lu(554)
         lu(560) = lu(560) - lu(468) * lu(554)
         lu(561) = lu(561) - lu(469) * lu(554)
         lu(562) = lu(562) - lu(470) * lu(554)
         lu(563) = lu(563) - lu(471) * lu(554)
         lu(564) = lu(564) - lu(472) * lu(554)
         lu(565) = lu(565) - lu(473) * lu(554)
         lu(630) = lu(630) - lu(465) * lu(629)
         lu(633) = lu(633) - lu(466) * lu(629)
         lu(634) = lu(634) - lu(467) * lu(629)
         lu(635) = lu(635) - lu(468) * lu(629)
         lu(636) = lu(636) - lu(469) * lu(629)
         lu(637) = lu(637) - lu(470) * lu(629)
         lu(638) = lu(638) - lu(471) * lu(629)
         lu(639) = lu(639) - lu(472) * lu(629)
         lu(640) = lu(640) - lu(473) * lu(629)
         lu(680) = lu(680) - lu(465) * lu(679)
         lu(683) = lu(683) - lu(466) * lu(679)
         lu(684) = lu(684) - lu(467) * lu(679)
         lu(685) = lu(685) - lu(468) * lu(679)
         lu(686) = lu(686) - lu(469) * lu(679)
         lu(687) = lu(687) - lu(470) * lu(679)
         lu(688) = lu(688) - lu(471) * lu(679)
         lu(689) = lu(689) - lu(472) * lu(679)
         lu(690) = lu(690) - lu(473) * lu(679)
         lu(710) = lu(710) - lu(465) * lu(709)
         lu(713) = lu(713) - lu(466) * lu(709)
         lu(714) = lu(714) - lu(467) * lu(709)
         lu(715) = lu(715) - lu(468) * lu(709)
         lu(716) = lu(716) - lu(469) * lu(709)
         lu(717) = lu(717) - lu(470) * lu(709)
         lu(718) = lu(718) - lu(471) * lu(709)
         lu(719) = lu(719) - lu(472) * lu(709)
         lu(720) = lu(720) - lu(473) * lu(709)
         lu(739) = lu(739) - lu(465) * lu(738)
         lu(742) = lu(742) - lu(466) * lu(738)
         lu(743) = lu(743) - lu(467) * lu(738)
         lu(744) = lu(744) - lu(468) * lu(738)
         lu(745) = lu(745) - lu(469) * lu(738)
         lu(746) = lu(746) - lu(470) * lu(738)
         lu(747) = lu(747) - lu(471) * lu(738)
         lu(748) = lu(748) - lu(472) * lu(738)
         lu(749) = lu(749) - lu(473) * lu(738)
         lu(792) = lu(792) - lu(465) * lu(791)
         lu(795) = lu(795) - lu(466) * lu(791)
         lu(796) = lu(796) - lu(467) * lu(791)
         lu(797) = lu(797) - lu(468) * lu(791)
         lu(798) = lu(798) - lu(469) * lu(791)
         lu(799) = lu(799) - lu(470) * lu(791)
         lu(800) = lu(800) - lu(471) * lu(791)
         lu(801) = lu(801) - lu(472) * lu(791)
         lu(802) = lu(802) - lu(473) * lu(791)
         lu(814) = lu(814) - lu(465) * lu(813)
         lu(817) = lu(817) - lu(466) * lu(813)
         lu(818) = lu(818) - lu(467) * lu(813)
         lu(819) = lu(819) - lu(468) * lu(813)
         lu(820) = lu(820) - lu(469) * lu(813)
         lu(821) = lu(821) - lu(470) * lu(813)
         lu(822) = lu(822) - lu(471) * lu(813)
         lu(823) = lu(823) - lu(472) * lu(813)
         lu(824) = lu(824) - lu(473) * lu(813)
         lu(482) = 1._r8 / lu(482)
         lu(483) = lu(483) * lu(482)
         lu(484) = lu(484) * lu(482)
         lu(485) = lu(485) * lu(482)
         lu(486) = lu(486) * lu(482)
         lu(487) = lu(487) * lu(482)
         lu(488) = lu(488) * lu(482)
         lu(489) = lu(489) * lu(482)
         lu(490) = lu(490) * lu(482)
         lu(502) = lu(502) - lu(483) * lu(500)
         lu(503) = lu(503) - lu(484) * lu(500)
         lu(504) = lu(504) - lu(485) * lu(500)
         lu(505) = lu(505) - lu(486) * lu(500)
         lu(507) = lu(507) - lu(487) * lu(500)
         lu(508) = lu(508) - lu(488) * lu(500)
         lu(509) = lu(509) - lu(489) * lu(500)
         lu(510) = lu(510) - lu(490) * lu(500)
         lu(557) = lu(557) - lu(483) * lu(555)
         lu(558) = lu(558) - lu(484) * lu(555)
         lu(559) = lu(559) - lu(485) * lu(555)
         lu(560) = lu(560) - lu(486) * lu(555)
         lu(562) = lu(562) - lu(487) * lu(555)
         lu(563) = lu(563) - lu(488) * lu(555)
         lu(564) = lu(564) - lu(489) * lu(555)
         lu(565) = lu(565) - lu(490) * lu(555)
         lu(632) = lu(632) - lu(483) * lu(630)
         lu(633) = lu(633) - lu(484) * lu(630)
         lu(634) = lu(634) - lu(485) * lu(630)
         lu(635) = lu(635) - lu(486) * lu(630)
         lu(637) = lu(637) - lu(487) * lu(630)
         lu(638) = lu(638) - lu(488) * lu(630)
         lu(639) = lu(639) - lu(489) * lu(630)
         lu(640) = lu(640) - lu(490) * lu(630)
         lu(682) = lu(682) - lu(483) * lu(680)
         lu(683) = lu(683) - lu(484) * lu(680)
         lu(684) = lu(684) - lu(485) * lu(680)
         lu(685) = lu(685) - lu(486) * lu(680)
         lu(687) = lu(687) - lu(487) * lu(680)
         lu(688) = lu(688) - lu(488) * lu(680)
         lu(689) = lu(689) - lu(489) * lu(680)
         lu(690) = lu(690) - lu(490) * lu(680)
         lu(712) = lu(712) - lu(483) * lu(710)
         lu(713) = lu(713) - lu(484) * lu(710)
         lu(714) = lu(714) - lu(485) * lu(710)
         lu(715) = lu(715) - lu(486) * lu(710)
         lu(717) = lu(717) - lu(487) * lu(710)
         lu(718) = lu(718) - lu(488) * lu(710)
         lu(719) = lu(719) - lu(489) * lu(710)
         lu(720) = lu(720) - lu(490) * lu(710)
         lu(741) = lu(741) - lu(483) * lu(739)
         lu(742) = lu(742) - lu(484) * lu(739)
         lu(743) = lu(743) - lu(485) * lu(739)
         lu(744) = lu(744) - lu(486) * lu(739)
         lu(746) = lu(746) - lu(487) * lu(739)
         lu(747) = lu(747) - lu(488) * lu(739)
         lu(748) = lu(748) - lu(489) * lu(739)
         lu(749) = lu(749) - lu(490) * lu(739)
         lu(794) = lu(794) - lu(483) * lu(792)
         lu(795) = lu(795) - lu(484) * lu(792)
         lu(796) = lu(796) - lu(485) * lu(792)
         lu(797) = lu(797) - lu(486) * lu(792)
         lu(799) = lu(799) - lu(487) * lu(792)
         lu(800) = lu(800) - lu(488) * lu(792)
         lu(801) = lu(801) - lu(489) * lu(792)
         lu(802) = lu(802) - lu(490) * lu(792)
         lu(816) = lu(816) - lu(483) * lu(814)
         lu(817) = lu(817) - lu(484) * lu(814)
         lu(818) = lu(818) - lu(485) * lu(814)
         lu(819) = lu(819) - lu(486) * lu(814)
         lu(821) = lu(821) - lu(487) * lu(814)
         lu(822) = lu(822) - lu(488) * lu(814)
         lu(823) = lu(823) - lu(489) * lu(814)
         lu(824) = lu(824) - lu(490) * lu(814)
         lu(501) = 1._r8 / lu(501)
         lu(502) = lu(502) * lu(501)
         lu(503) = lu(503) * lu(501)
         lu(504) = lu(504) * lu(501)
         lu(505) = lu(505) * lu(501)
         lu(506) = lu(506) * lu(501)
         lu(507) = lu(507) * lu(501)
         lu(508) = lu(508) * lu(501)
         lu(509) = lu(509) * lu(501)
         lu(510) = lu(510) * lu(501)
         lu(557) = lu(557) - lu(502) * lu(556)
         lu(558) = lu(558) - lu(503) * lu(556)
         lu(559) = lu(559) - lu(504) * lu(556)
         lu(560) = lu(560) - lu(505) * lu(556)
         lu(561) = lu(561) - lu(506) * lu(556)
         lu(562) = lu(562) - lu(507) * lu(556)
         lu(563) = lu(563) - lu(508) * lu(556)
         lu(564) = lu(564) - lu(509) * lu(556)
         lu(565) = lu(565) - lu(510) * lu(556)
         lu(632) = lu(632) - lu(502) * lu(631)
         lu(633) = lu(633) - lu(503) * lu(631)
         lu(634) = lu(634) - lu(504) * lu(631)
         lu(635) = lu(635) - lu(505) * lu(631)
         lu(636) = lu(636) - lu(506) * lu(631)
         lu(637) = lu(637) - lu(507) * lu(631)
         lu(638) = lu(638) - lu(508) * lu(631)
         lu(639) = lu(639) - lu(509) * lu(631)
         lu(640) = lu(640) - lu(510) * lu(631)
         lu(682) = lu(682) - lu(502) * lu(681)
         lu(683) = lu(683) - lu(503) * lu(681)
         lu(684) = lu(684) - lu(504) * lu(681)
         lu(685) = lu(685) - lu(505) * lu(681)
         lu(686) = lu(686) - lu(506) * lu(681)
         lu(687) = lu(687) - lu(507) * lu(681)
         lu(688) = lu(688) - lu(508) * lu(681)
         lu(689) = lu(689) - lu(509) * lu(681)
         lu(690) = lu(690) - lu(510) * lu(681)
         lu(712) = lu(712) - lu(502) * lu(711)
         lu(713) = lu(713) - lu(503) * lu(711)
         lu(714) = lu(714) - lu(504) * lu(711)
         lu(715) = lu(715) - lu(505) * lu(711)
         lu(716) = lu(716) - lu(506) * lu(711)
         lu(717) = lu(717) - lu(507) * lu(711)
         lu(718) = lu(718) - lu(508) * lu(711)
         lu(719) = lu(719) - lu(509) * lu(711)
         lu(720) = lu(720) - lu(510) * lu(711)
         lu(741) = lu(741) - lu(502) * lu(740)
         lu(742) = lu(742) - lu(503) * lu(740)
         lu(743) = lu(743) - lu(504) * lu(740)
         lu(744) = lu(744) - lu(505) * lu(740)
         lu(745) = lu(745) - lu(506) * lu(740)
         lu(746) = lu(746) - lu(507) * lu(740)
         lu(747) = lu(747) - lu(508) * lu(740)
         lu(748) = lu(748) - lu(509) * lu(740)
         lu(749) = lu(749) - lu(510) * lu(740)
         lu(794) = lu(794) - lu(502) * lu(793)
         lu(795) = lu(795) - lu(503) * lu(793)
         lu(796) = lu(796) - lu(504) * lu(793)
         lu(797) = lu(797) - lu(505) * lu(793)
         lu(798) = lu(798) - lu(506) * lu(793)
         lu(799) = lu(799) - lu(507) * lu(793)
         lu(800) = lu(800) - lu(508) * lu(793)
         lu(801) = lu(801) - lu(509) * lu(793)
         lu(802) = lu(802) - lu(510) * lu(793)
         lu(816) = lu(816) - lu(502) * lu(815)
         lu(817) = lu(817) - lu(503) * lu(815)
         lu(818) = lu(818) - lu(504) * lu(815)
         lu(819) = lu(819) - lu(505) * lu(815)
         lu(820) = lu(820) - lu(506) * lu(815)
         lu(821) = lu(821) - lu(507) * lu(815)
         lu(822) = lu(822) - lu(508) * lu(815)
         lu(823) = lu(823) - lu(509) * lu(815)
         lu(824) = lu(824) - lu(510) * lu(815)
         lu(517) = 1._r8 / lu(517)
         lu(518) = lu(518) * lu(517)
         lu(519) = lu(519) * lu(517)
         lu(520) = lu(520) * lu(517)
         lu(521) = lu(521) * lu(517)
         lu(522) = lu(522) * lu(517)
         lu(523) = lu(523) * lu(517)
         lu(524) = lu(524) * lu(517)
         lu(525) = lu(525) * lu(517)
         lu(558) = lu(558) - lu(518) * lu(557)
         lu(559) = lu(559) - lu(519) * lu(557)
         lu(560) = lu(560) - lu(520) * lu(557)
         lu(561) = lu(561) - lu(521) * lu(557)
         lu(562) = lu(562) - lu(522) * lu(557)
         lu(563) = lu(563) - lu(523) * lu(557)
         lu(564) = lu(564) - lu(524) * lu(557)
         lu(565) = lu(565) - lu(525) * lu(557)
         lu(633) = lu(633) - lu(518) * lu(632)
         lu(634) = lu(634) - lu(519) * lu(632)
         lu(635) = lu(635) - lu(520) * lu(632)
         lu(636) = lu(636) - lu(521) * lu(632)
         lu(637) = lu(637) - lu(522) * lu(632)
         lu(638) = lu(638) - lu(523) * lu(632)
         lu(639) = lu(639) - lu(524) * lu(632)
         lu(640) = lu(640) - lu(525) * lu(632)
         lu(683) = lu(683) - lu(518) * lu(682)
         lu(684) = lu(684) - lu(519) * lu(682)
         lu(685) = lu(685) - lu(520) * lu(682)
         lu(686) = lu(686) - lu(521) * lu(682)
         lu(687) = lu(687) - lu(522) * lu(682)
         lu(688) = lu(688) - lu(523) * lu(682)
         lu(689) = lu(689) - lu(524) * lu(682)
         lu(690) = lu(690) - lu(525) * lu(682)
         lu(713) = lu(713) - lu(518) * lu(712)
         lu(714) = lu(714) - lu(519) * lu(712)
         lu(715) = lu(715) - lu(520) * lu(712)
         lu(716) = lu(716) - lu(521) * lu(712)
         lu(717) = lu(717) - lu(522) * lu(712)
         lu(718) = lu(718) - lu(523) * lu(712)
         lu(719) = lu(719) - lu(524) * lu(712)
         lu(720) = lu(720) - lu(525) * lu(712)
         lu(742) = lu(742) - lu(518) * lu(741)
         lu(743) = lu(743) - lu(519) * lu(741)
         lu(744) = lu(744) - lu(520) * lu(741)
         lu(745) = lu(745) - lu(521) * lu(741)
         lu(746) = lu(746) - lu(522) * lu(741)
         lu(747) = lu(747) - lu(523) * lu(741)
         lu(748) = lu(748) - lu(524) * lu(741)
         lu(749) = lu(749) - lu(525) * lu(741)
         lu(762) = lu(762) - lu(518) * lu(761)
         lu(763) = lu(763) - lu(519) * lu(761)
         lu(764) = lu(764) - lu(520) * lu(761)
         lu(765) = lu(765) - lu(521) * lu(761)
         lu(766) = lu(766) - lu(522) * lu(761)
         lu(767) = lu(767) - lu(523) * lu(761)
         lu(768) = lu(768) - lu(524) * lu(761)
         lu(769) = lu(769) - lu(525) * lu(761)
         lu(795) = lu(795) - lu(518) * lu(794)
         lu(796) = lu(796) - lu(519) * lu(794)
         lu(797) = lu(797) - lu(520) * lu(794)
         lu(798) = lu(798) - lu(521) * lu(794)
         lu(799) = lu(799) - lu(522) * lu(794)
         lu(800) = lu(800) - lu(523) * lu(794)
         lu(801) = lu(801) - lu(524) * lu(794)
         lu(802) = lu(802) - lu(525) * lu(794)
         lu(817) = lu(817) - lu(518) * lu(816)
         lu(818) = lu(818) - lu(519) * lu(816)
         lu(819) = lu(819) - lu(520) * lu(816)
         lu(820) = lu(820) - lu(521) * lu(816)
         lu(821) = lu(821) - lu(522) * lu(816)
         lu(822) = lu(822) - lu(523) * lu(816)
         lu(823) = lu(823) - lu(524) * lu(816)
         lu(824) = lu(824) - lu(525) * lu(816)
         lu(558) = 1._r8 / lu(558)
         lu(559) = lu(559) * lu(558)
         lu(560) = lu(560) * lu(558)
         lu(561) = lu(561) * lu(558)
         lu(562) = lu(562) * lu(558)
         lu(563) = lu(563) * lu(558)
         lu(564) = lu(564) * lu(558)
         lu(565) = lu(565) * lu(558)
         lu(634) = lu(634) - lu(559) * lu(633)
         lu(635) = lu(635) - lu(560) * lu(633)
         lu(636) = lu(636) - lu(561) * lu(633)
         lu(637) = lu(637) - lu(562) * lu(633)
         lu(638) = lu(638) - lu(563) * lu(633)
         lu(639) = lu(639) - lu(564) * lu(633)
         lu(640) = lu(640) - lu(565) * lu(633)
         lu(684) = lu(684) - lu(559) * lu(683)
         lu(685) = lu(685) - lu(560) * lu(683)
         lu(686) = lu(686) - lu(561) * lu(683)
         lu(687) = lu(687) - lu(562) * lu(683)
         lu(688) = lu(688) - lu(563) * lu(683)
         lu(689) = lu(689) - lu(564) * lu(683)
         lu(690) = lu(690) - lu(565) * lu(683)
         lu(714) = lu(714) - lu(559) * lu(713)
         lu(715) = lu(715) - lu(560) * lu(713)
         lu(716) = lu(716) - lu(561) * lu(713)
         lu(717) = lu(717) - lu(562) * lu(713)
         lu(718) = lu(718) - lu(563) * lu(713)
         lu(719) = lu(719) - lu(564) * lu(713)
         lu(720) = lu(720) - lu(565) * lu(713)
         lu(743) = lu(743) - lu(559) * lu(742)
         lu(744) = lu(744) - lu(560) * lu(742)
         lu(745) = lu(745) - lu(561) * lu(742)
         lu(746) = lu(746) - lu(562) * lu(742)
         lu(747) = lu(747) - lu(563) * lu(742)
         lu(748) = lu(748) - lu(564) * lu(742)
         lu(749) = lu(749) - lu(565) * lu(742)
         lu(763) = lu(763) - lu(559) * lu(762)
         lu(764) = lu(764) - lu(560) * lu(762)
         lu(765) = lu(765) - lu(561) * lu(762)
         lu(766) = lu(766) - lu(562) * lu(762)
         lu(767) = lu(767) - lu(563) * lu(762)
         lu(768) = lu(768) - lu(564) * lu(762)
         lu(769) = lu(769) - lu(565) * lu(762)
         lu(796) = lu(796) - lu(559) * lu(795)
         lu(797) = lu(797) - lu(560) * lu(795)
         lu(798) = lu(798) - lu(561) * lu(795)
         lu(799) = lu(799) - lu(562) * lu(795)
         lu(800) = lu(800) - lu(563) * lu(795)
         lu(801) = lu(801) - lu(564) * lu(795)
         lu(802) = lu(802) - lu(565) * lu(795)
         lu(818) = lu(818) - lu(559) * lu(817)
         lu(819) = lu(819) - lu(560) * lu(817)
         lu(820) = lu(820) - lu(561) * lu(817)
         lu(821) = lu(821) - lu(562) * lu(817)
         lu(822) = lu(822) - lu(563) * lu(817)
         lu(823) = lu(823) - lu(564) * lu(817)
         lu(824) = lu(824) - lu(565) * lu(817)
      end subroutine lu_fac10
      subroutine lu_fac11( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
         lu(634) = 1._r8 / lu(634)
         lu(635) = lu(635) * lu(634)
         lu(636) = lu(636) * lu(634)
         lu(637) = lu(637) * lu(634)
         lu(638) = lu(638) * lu(634)
         lu(639) = lu(639) * lu(634)
         lu(640) = lu(640) * lu(634)
         lu(685) = lu(685) - lu(635) * lu(684)
         lu(686) = lu(686) - lu(636) * lu(684)
         lu(687) = lu(687) - lu(637) * lu(684)
         lu(688) = lu(688) - lu(638) * lu(684)
         lu(689) = lu(689) - lu(639) * lu(684)
         lu(690) = lu(690) - lu(640) * lu(684)
         lu(715) = lu(715) - lu(635) * lu(714)
         lu(716) = lu(716) - lu(636) * lu(714)
         lu(717) = lu(717) - lu(637) * lu(714)
         lu(718) = lu(718) - lu(638) * lu(714)
         lu(719) = lu(719) - lu(639) * lu(714)
         lu(720) = lu(720) - lu(640) * lu(714)
         lu(744) = lu(744) - lu(635) * lu(743)
         lu(745) = lu(745) - lu(636) * lu(743)
         lu(746) = lu(746) - lu(637) * lu(743)
         lu(747) = lu(747) - lu(638) * lu(743)
         lu(748) = lu(748) - lu(639) * lu(743)
         lu(749) = lu(749) - lu(640) * lu(743)
         lu(764) = lu(764) - lu(635) * lu(763)
         lu(765) = lu(765) - lu(636) * lu(763)
         lu(766) = lu(766) - lu(637) * lu(763)
         lu(767) = lu(767) - lu(638) * lu(763)
         lu(768) = lu(768) - lu(639) * lu(763)
         lu(769) = lu(769) - lu(640) * lu(763)
         lu(797) = lu(797) - lu(635) * lu(796)
         lu(798) = lu(798) - lu(636) * lu(796)
         lu(799) = lu(799) - lu(637) * lu(796)
         lu(800) = lu(800) - lu(638) * lu(796)
         lu(801) = lu(801) - lu(639) * lu(796)
         lu(802) = lu(802) - lu(640) * lu(796)
         lu(819) = lu(819) - lu(635) * lu(818)
         lu(820) = lu(820) - lu(636) * lu(818)
         lu(821) = lu(821) - lu(637) * lu(818)
         lu(822) = lu(822) - lu(638) * lu(818)
         lu(823) = lu(823) - lu(639) * lu(818)
         lu(824) = lu(824) - lu(640) * lu(818)
         lu(685) = 1._r8 / lu(685)
         lu(686) = lu(686) * lu(685)
         lu(687) = lu(687) * lu(685)
         lu(688) = lu(688) * lu(685)
         lu(689) = lu(689) * lu(685)
         lu(690) = lu(690) * lu(685)
         lu(716) = lu(716) - lu(686) * lu(715)
         lu(717) = lu(717) - lu(687) * lu(715)
         lu(718) = lu(718) - lu(688) * lu(715)
         lu(719) = lu(719) - lu(689) * lu(715)
         lu(720) = lu(720) - lu(690) * lu(715)
         lu(745) = lu(745) - lu(686) * lu(744)
         lu(746) = lu(746) - lu(687) * lu(744)
         lu(747) = lu(747) - lu(688) * lu(744)
         lu(748) = lu(748) - lu(689) * lu(744)
         lu(749) = lu(749) - lu(690) * lu(744)
         lu(765) = lu(765) - lu(686) * lu(764)
         lu(766) = lu(766) - lu(687) * lu(764)
         lu(767) = lu(767) - lu(688) * lu(764)
         lu(768) = lu(768) - lu(689) * lu(764)
         lu(769) = lu(769) - lu(690) * lu(764)
         lu(798) = lu(798) - lu(686) * lu(797)
         lu(799) = lu(799) - lu(687) * lu(797)
         lu(800) = lu(800) - lu(688) * lu(797)
         lu(801) = lu(801) - lu(689) * lu(797)
         lu(802) = lu(802) - lu(690) * lu(797)
         lu(820) = lu(820) - lu(686) * lu(819)
         lu(821) = lu(821) - lu(687) * lu(819)
         lu(822) = lu(822) - lu(688) * lu(819)
         lu(823) = lu(823) - lu(689) * lu(819)
         lu(824) = lu(824) - lu(690) * lu(819)
         lu(716) = 1._r8 / lu(716)
         lu(717) = lu(717) * lu(716)
         lu(718) = lu(718) * lu(716)
         lu(719) = lu(719) * lu(716)
         lu(720) = lu(720) * lu(716)
         lu(746) = lu(746) - lu(717) * lu(745)
         lu(747) = lu(747) - lu(718) * lu(745)
         lu(748) = lu(748) - lu(719) * lu(745)
         lu(749) = lu(749) - lu(720) * lu(745)
         lu(766) = lu(766) - lu(717) * lu(765)
         lu(767) = lu(767) - lu(718) * lu(765)
         lu(768) = lu(768) - lu(719) * lu(765)
         lu(769) = lu(769) - lu(720) * lu(765)
         lu(799) = lu(799) - lu(717) * lu(798)
         lu(800) = lu(800) - lu(718) * lu(798)
         lu(801) = lu(801) - lu(719) * lu(798)
         lu(802) = lu(802) - lu(720) * lu(798)
         lu(821) = lu(821) - lu(717) * lu(820)
         lu(822) = lu(822) - lu(718) * lu(820)
         lu(823) = lu(823) - lu(719) * lu(820)
         lu(824) = lu(824) - lu(720) * lu(820)
         lu(746) = 1._r8 / lu(746)
         lu(747) = lu(747) * lu(746)
         lu(748) = lu(748) * lu(746)
         lu(749) = lu(749) * lu(746)
         lu(767) = lu(767) - lu(747) * lu(766)
         lu(768) = lu(768) - lu(748) * lu(766)
         lu(769) = lu(769) - lu(749) * lu(766)
         lu(800) = lu(800) - lu(747) * lu(799)
         lu(801) = lu(801) - lu(748) * lu(799)
         lu(802) = lu(802) - lu(749) * lu(799)
         lu(822) = lu(822) - lu(747) * lu(821)
         lu(823) = lu(823) - lu(748) * lu(821)
         lu(824) = lu(824) - lu(749) * lu(821)
         lu(767) = 1._r8 / lu(767)
         lu(768) = lu(768) * lu(767)
         lu(769) = lu(769) * lu(767)
         lu(801) = lu(801) - lu(768) * lu(800)
         lu(802) = lu(802) - lu(769) * lu(800)
         lu(823) = lu(823) - lu(768) * lu(822)
         lu(824) = lu(824) - lu(769) * lu(822)
         lu(801) = 1._r8 / lu(801)
         lu(802) = lu(802) * lu(801)
         lu(824) = lu(824) - lu(802) * lu(823)
         lu(824) = 1._r8 / lu(824)
      end subroutine lu_fac11
      subroutine lu_fac( lu )
      use shr_kind_mod, only : r8 => shr_kind_r8
      implicit none
!-----------------------------------------------------------------------
! ... dummy args
!-----------------------------------------------------------------------
      real(r8), intent(inout) :: lu(:)
      call lu_fac01( lu )
      call lu_fac02( lu )
      call lu_fac03( lu )
      call lu_fac04( lu )
      call lu_fac05( lu )
      call lu_fac06( lu )
      call lu_fac07( lu )
      call lu_fac08( lu )
      call lu_fac09( lu )
      call lu_fac10( lu )
      call lu_fac11( lu )
      end subroutine lu_fac
      end module mo_lu_factor
